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The selfconsistent cranking approach is extended to the 
case of rotation about an axis which is tilted with respect 
to the principal axes of the deformed potential (Tilted Axis 
Cranking). Expressions for the energies and the intra bands 
electromagnetic transition probabilities are given. The mean 
field solutions are interpreted in terms of quantal rotational 
states. The construction of the quasi particle configurations 
and the elimination of spurious states is discussed. The ap- 
plication of the theory to high spin data is demonstrated by 
analyzing the multi quasi particle bands in the nuclide-s with 
N = 102, 103 and Z = 71, 72, 73. 

PACS numbers: 21.60.-ii 

Keywords: High- if rotational bands, Tilted Axis Crank- 
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I. INTRODUCTION 

Since its introduction in ref . [Q , the Tilted Axis Crank- 
ing (TAC) approach has turned out to be quite successful 
in describing AI — 1 rotational bands In partic- 

ular it has led to the understanding of the appearance 
of regular magnetic dipole bands in nearly spherical nu- 
clei [^-^ . The physical aspects of these investigations 
have been reviewed in Though different aspects 

of the actual calculations were touched in these papers, 
a comprehensive presentation of the theory, the calcu- 
lational methods, and the practical application of the 
TAC approach is still missing. In the present paper we 
provide it by using the rotational bands in the nuclides 
with N = 102, 103 and Z = 71, 72, 73 as examples. The 
TAC approach has been applied so far only for potentials 
of the Nilsson type, which are combined with a pairing 
plus quadrupole model interaction or the shell correction 
method for finding the deformation. A systematic expo- 
sure of the applied techniques, the experiences gathered 
as well as the successes and limitations of the TAC ap- 
proach as it stands seems to be timely for two reasons. 
On the one hand hand it is meant as guideline for appli- 
cation of the existing program system, which has turned 
out quite useful in the data analysis. One the other hand 
it may serve as a starting point for more sophisticated 
versions of the TAC mean field approximation, as the up 
to date versions of the Hartree-Fock approximation or 
the Relativistic Mean Field approach. 

The earliest invocations of cranking about a non- 
principal axis were in the context of wobbling motion 
p7|-p0[. Kerman and Onishi |Q pointed out the pos- 
sibility of uniform rotation about a non-principal axis. 



Frisk and Bengtsson |31[] demonstrated the existence of 
such solutions for realistic nuclei and discussed conditions 
where to expect them |^,^. Goodman jsj] demon- 
strated that the moments of inertia may strongly depend 
on the orientation of the rotational axis, which implies 
the possibility of uniform rotation about a tilted axis. 
However, these studies did not give the physical interpre- 
tation of the TAC solutions and left open the question if 
taking into account the self-consistency with respect to 
the shape degrees of freedom would not result in rotation 
about a principal axis. In fact, the investigation of the 
rotating harmonic oscillator by Cuypers |35| and a few 
level model by Nazarewicz and Szymanski |M| seemed to 
support the latter possibility. Frauendorf ^ found the 
first fully self-consistent TAC solutions and gave their in- 
terpretation in terms of AI = 1 rotational bands. This 
marked the origin of the fully fledged Tilted Axis Crank- 
ing (TAC) ap pro ach. 

Marshalek [ p7|j3q| studied the possibility of tilted rota- 
tion generated by superpositions of collective vibrations, 
while Alhassid and Bush , Goodman , and Dodaro 
and Goodman included the tilt of the rotational 

axis in their analysis of nuclei at nonzero temperature. 
A recent reinvestigation of the rotating harmonic oscilla- 
tor by Heiss and Nazmitdinov jij] claims the existence 
of TAC solutions within this model, in contrast to [ p5[ . 
Horibata and Onishi Horibata et al. and Donau 
et al. |4^] have started to investigate the dynamics of the 
orientation angles in the frame of the Generator Coordi- 
nate Method. 

Section ^ presents the relevant expressions for the en- 
ergies and electro-magnetic transition matrix elements. 
It discusses the interpretation of the cranking solutions, 
important technical aspects and approximations that 
help to find the solution of the TAC equations in an ef- 
ficient way. It investigates the relation of TAC to the 
treatment of AI = 1 bands in the framework of the stan- 
dard Principal Axis Cranking (PAC) approach. It ex- 
plains how to read the quasi particle diagrams. Section 
III analyses the rotational bands in the yrast region of the 
nuchdes with N = 102, 103 and Z = 71, 72, 73. The main 
purpose is to illustrate how to construct the multi quasi 
particle configurations and how to relate them to the ex- 
perimental rotational bands. Merits and limitations of 
the method will be exposed and compared with the stan- 
dard Cranking approach. We are not going to optimize 
all parameters of the mean field for each configuration. In 
the spirit of the Cranked Shell Model approach only 
semi quantitative agreement with the data is sought, the 
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focus being the quahtative structure of the band spec- 
trum. Well deformed nuclei are taken as examples be- 
cause the assumption of one and the same deformation 
for the various quasi particle configurations is realistic. 
The specific nuclcon numbers are chosen because a large 
number of high K bands and low K bands have been 
found in these nuclides. This makes them an appropriate 
test ground for the TAG model. This paper is restricted 
to the HFB approximation for pairing. A more sophis- 
ticated version of TAG based on particle number projec- 
tion will be presented separately |^^. Since the change 
of the pair field is not in the concern of this paper but 
rather an unwanted complication, the self-consistency of 
this degree of freedom is treated in a schematic way. 

Section |^ provides a set of rules for the analysis of 
rotational bands in terms TAG, which is meant as a ref- 
erence for potential users of the method. The conclusion 
are given in section 



II. TILTED AXIS CRANKING 

A. General layout 

Two versions of the TAG have been developed and ap- 
plied 

i) The pairing plus quadrupole model (PQTAG) ^ 

ii) Shell correction method (SGTAG) 

In the subsecti ons |IIB| - p^Il| we present the PQTAG 
in detail. Section |[IE| describes the differences between 
SGTAG and PQTAG. The PQTAG it more appropri- 
ate for small deformations, whereas SG TAG is better 
suited for large deformation. Subsection [IJ discusses 



the schematic treatment of pairing and subsection II K 
explains the how to read the quasi particle diagrams. 



B. The pairing plus quadrupole model (PQTAG) 

We assume that the rotational axis is the z-axis and 
start with the two-body Routhian 



H' = H ~ ujJ^ 



(1) 



It consists of the rotationally invariant two-body Hamil- 
tonian H and the constraint uoJz which ensures that the 
low-lying states have a finite angular momentum projec- 
tion J =< J2 >. As a two body Hamiltonian, the pairing 
plus quadrupole interaction is used. 



H ^ H. 



sph 



XN. 



(2) 



The model and its justification are described in the text- 
books (see, for example. Ring and Schuck ^]). We use a 
slightly modified version, which is constructed such that 
the derived mean field Hamiltonian coincides with the 



popular Nilsson Hamiltonian (see, for example. Ring and 
Schuck and Nilsson and Ragnarsson ||5^). The mo- 
tivation is that the parameters of the Nilsson Hamilto- 
nian are carefully adjusted and that it is useful to have 
the same standard mean field for nuclei with large defor- 
mation, where the shell correction method pl| is more 
appropriate. Thus, the spherical part 



■ph 



k 



(3) 



is parameterized in the same way as the Nilsson Hamil- 
tonian. For the calculations in this paper we use the set 
of parameters given in . 

The pairing interaction is defined by the monopole pair 
field 



E 

fc>0 



A^^k^^k- 



(4) 



Here k is the time reversed state of k. The quadrupole 
interaction is defined by the operators [] 



/ 4-7r 

= E V — < ^I'^'^aJfc' > 



(5) 



fefe' 



In order to simplify the notation all expressions are writ- 
ten only for one kind of particles. They are understood 
as sums of a proton and a neutron part. 

The wave function is approximated by the Hartree - 
Fock - Bogoljubov (HFB) mean field expression | >. Ne- 
glecting exchange terms, the HFB - Routhian becomes 



- UjJz 



(6) 



The self-consistency equations determine the deformed 
part of the potential 



9m = X < > 
and the pair potential 

A = G < P > . 



(7) 



(8) 



The chemical potential A is fixed by the standard condi- 
tion 



N ^< N > 
The quasi particle operators 



(9) 



^This definition of the quadrupole operators corresponds to 
Qo = r^P2(cosi9), with P2 being the Legendre polynomial. 
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obey the equations of motion 
[h', a 



(10) 



(11) 



which define the well known HFB eigenvalue equations 
for the quasi particle amplitudes Uki and Vki- The ex- 
plicit form of these equations can be found in . They 
have the familiar symmetry under particle hole conju- 
gation, which has the consequence that for each quasi 
particle solution i there is a conjugate z+ with 



Uki+ = Vfej, Vki+ = Uk 



(12) 



The quasi particles have good parity, but in general no 
good signature. The consequences of good or broken sig- 
nature will be discussed in subsection II C . 



The quasi particle operators refer to the vacuum state 
|0 >, which is defined by the condition 

a, |0>=OVi. (13) 

They define the excited quasi particle configurations 



12, ■■■ >^aj^aj_^ 



10 > 



(14) 



The rules and strategies for constructing quasi particle 
configurations from them will be be discussed below by 
means of concrete examples. 

The set of HFB eq. (|) - (|l|) can be solved for any con- 
figuration 1*1, 12 ■ ■ ■ >■ For the self-consistent solution, 
the total Routhian 



E' =< H' > 



has an extremum 



9m 



(15) 



(16) 



The total energy as function of the angular momentum 
J J. > (17) 

is given by 



E{J) = E'{lj) +ujJ{lu) 



(18) 



where the eq. ( p7|) implicitly fixes w( J). The total energy 
is also extremal for the self-consistent solution 



^^1 -0 ^^1 -0 

^U- 0,-1,-0, 



(19) 



where the derivatives have to be taken at a fixed value of 
J. For a family of self-consistent solutions \uj >, found for 
different values of oj, there hold the canonical relations 



dE' _ J dE _ ^ 
du> ' dJ 



(20) 



In ref. it was shown that for a self-consistent solu- 
tion the vector of the angular velocity 



LU = {uJx,UJy,UJz) = (0,0, w) 



(21) 



and the vector of the expectation values of the angular 
momentum components 



J =((4), (J,), (J.)) 



(22) 



must be parallel. The argument is as follows. Since the 
interaction is rotational invariant, one has 



(23) 
(24) 



{[H',Jy])^-tuj{J^). 

Since the left-hand sides are small variations of £", the 
stationarity of E' implies 



< > = < Jy >= 0, 



(25) 



1. e. 



uj\\J. (26) 
This holds also in the intrinsic frame of reference, which 



will be discussed in sect. II D 



C. Tilted Solutions 

The formalism presented above is the well known 
"Cranking Model" as laid out in the textbooks ||9|,|o|. 
The "Tilted- Axis Cranking" (TAC) version |] accounts 
for the possibility that the principal axes (PA) of the 
quadrupole tensor need not to coincide with the ro- 
tational axis (z). Hence, one has to distinguish between 
two possibilities: 

• Principal Axis Cranking (PAC). The z - axis (ro- 
tational or cranking axis) coincides with one of the 
PA. Then, the signature r is a good quantum num- 
ber, i. e. 



-ittJz 



TT, a, w >= r|7r, a, w > 



(27) 



Following we indicate the signature r — e~"" 
by the signature exponent a. The quasi particle 
configuration |7r, Q;,a; > describes a AJ = 2 rota- 
tional band, the spins of which take the values ]47| ] 



I = a + even number. 



(28) 



• Tilted Axis Cranking (TAC). The z - axis does not 
coincides with one of the PA, i. e. it is tilted away 
from the PA. Then, 



^ Ivr, u! e 



\ttuj > 



(29) 



The signature is no longer a good quantum number. 
The quasi particle configuration |7r, to > describes 
a A/ = 1 rotational band of given parity. 
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The different interpretation of solutions with different 
symmetry is characteristic for spontaneous symmetry 
breaking in the mean field approximation. It makes it 
necessary to eliminate spurious configurations and will 
lead to discontinuities when the symmetry changes as a 
function of the frequency uj. These problems which will 
be discussed in the subsections [II B by means of concrete 
examples. 



D. Intrinsic coordinates 

It is useful to reformulate the TAG approach in the 
frame of the PA of the quadrupole tensor. This "intrin- 
sic" coordinate system is defined by demanding that the 
components of the quadrupole tensor satisfy the condi- 
tions 



9-1 9i = 0, g_2 



92- 



(30) 



The orientation of the PA, which are denoted by 1, 2, and 
3, with respect to the lab frame is fixed by the three Euler 
angles ip^ d and (f, the meaning of which is illustrated in 
fig. |l|. The two intrinsic quadrupole moments and q'2 
specify the deformation of the potential. The quadrupole 
moments in the lab frame are related to them by 



D 



(31) 



where -D^^(?A, ip) are the Wigner _D-functions ^. 

The different angles correspond to one and the same 
intrinsic state. They are degenerate. We choose the one 
with ^ = 0. We restrict the consideration to planar 
TAG solutions, which is the case when the z-axis lies in 
one of the three principal planes defined by the PA. We 
assume that it lies in the 1-3 plane, i. e. choose f — Q. 
In the case of axial shapes, this is one choice from the 
equivalent solutions differing by the angle Lp. For triaxial 
shapes one may always relabel the PA by means of the 
shape parameterization, letting the triaxiality parameter 
7 vary within an interval of 180° |50|. Which axes of 



the triaxial potential lie in the 1-3 plane can be found in 
table |. It is seen that all three possibilities appear in 
the half-plane -60° < 7 < 120°. The other half-plane is 
a repetition with the axes 1 and 3 exchanged. 

With the above mentioned restrictions and conven- 
tions the deformed potential is fixed by the two intrin- 
sic quadrupole moments q^ and q2 and the orientation 
("tilt") angle d between the 3 - and the z - axis, which is 
the direction of the rotational axis. In the intrinsic frame 
the HFB Routhian reads 



h' = KpH - q'oQ'o - 92(^2 + QI2) - A(P+ + P) 

-AA^-w(sin z9 Ji -l-cos 1? J3). (32) 

Figs. H and ^ show examples of the the quasi particle 
levels as functions of the rotational frequency lu and 
the orientation angle 1). 

The shape is fixed by the two equations 



q'o ^ K < Q'q >, q^ ^ K < Q2 > 



(33) 



and the orientation angle ?? by the condition that the 
expectation value of the angular momentum and the an- 
gular velocity must have the same direction, i. e. 

J= (< Ji >, 0, < J3 >) II uj ^ (w sin -d, 0, uj cos d), 

(34) 

respectively. These parameters correspond to extrema of 
total Routhian, that is 



= 0, 



= 0, 



I2 



dE' 



= 0. 



(35) 



Of course only the minima are interpreted as bands. 

In praxis it is convenient to solve the equation (jsj) 
for each combination of q^ and ^2 which is needed to 
obtain the shape from the equations (^) with the desired 
accuracy. Very often it is enough to determine the shape 
for one value of uj and then keep it fixed for other values, 
only calculating the orientation angle by means of the 
condition (|^. 

Using the Gartesian representation of the quadrupole 
moments, the HFB Routhian (21) becomes the modified 
oscillator potential |[l9l|5C|] 



iy=l 
5+ 



2M 



-A(P+ + P) - A7V - uj{sin i)Ji+ cos § J3), (36) 

where the oscillator frequencies are parameterized by 
means of Nilsson's deformation parameters 5 and 7 (cf. 
e. g. |0|), 



^00 



1 2. , 2^v- 
1 - -6 cos(7 - — ) 



(37) 



The only difference to the standard modified oscillator 
model is that there is no volume conservation in the pair- 
ing plus quadrupole model. Since we are only interested 
in small deformation the coupling between the oscillator 
shells is not taken into account when diagonalizing the 
HFB Routhian (^6|). Solving the self-consistency equa- 
tion (|33| ) and calculating the total energy, the coupling 
between the oscillator shells is also neglected. 



The convention of 



is used. 
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E. Strutinsky Renormalization (SCTAC) 



F. Electro - magnetic matrix elements 



An alternative version of TAG starts with the modi- 
fied oscihator Routhian (^. As e. g. described in Q, 
stretched coordinates are introduced and the matrix el- 
ements < -/V| J^I-ZV ± 2 > are neglected in the stretched 
basis. This is a standard procedure which takes into ac- 
count most of the couplings between the oscillator shells. 
The oscillator frequencies are parameterized by means of 
Nilsson's alternative set of deformation parameters e and 
7 , 



1 - -e cos(7 - — ) 



(38) 



where the condition of volume conservation cjq — (^1(^2^3 
fixes . The total Routhian is obtained by applying 
the Strutinsky renormalization to the energy of the non- 
rotating system Eq. This kind of approach has turned 
out to be a quite reliable calculation scheme in the case 
of standard PAC . One minimizes the total Routhian^ 

E'{uj, d, e, £4, 7, A, A) = Eld{£, £4, 7) - E{e, £4, 7) 

+ <h' > +(2A - G < P >) < P >, (39) 

where | >= |w, £, £4, 7, A, A > is a quasi parti- 
cle configuration belonging to the mean field Routhian 
/i'([j, t^, £, £4, 7, A, A) as defined above. The smooth en- 
ergy E is calculated from the single-particle energies, 
which are the eigenvalues of h'{uj = 0, = 0, £, £4, 7, A = 
0) by means of Strutinsky averaging The expres- 

sions for the liquid drop energy _Ei£)(£, £4, 7) are given, 
for example, in |^o| , where also the averaging procedure 
is described. For given £, £4 and 7, the tilt angle is 
determined by means of the condition (0). Then, the 
minimum of i?'(ti;, £, £4, 7) with respect to the deforma- 
tion parameters is found. Since \uj > is an eigenfunction 
of h'{uj, the Routhian < uj, "dlh'loj, i? > is stationary at 
the angle where the condition (|3J) is fulfilled and so is 
E' because the other terms do not depend on d. Hence, 
the procedure determines a stationary point with respect 
to the mean field parameters and the canonical relations 
( pO| ) are satisfied. 

The SCTAC approach is preferred to the PQTAC ver- 
sion for well deformed nuclei, because it is a reliable stan- 
dard method for determining large deformations. In the 
calculations of well deformed nuclei it is usually a good 
approximation to keep the deformations fixed within a 
rotational band. However this is a matter of the needed 
accuracy and of how much effort one is willing to invest. 



The intra band Ml - transition matrix element is cal- 
culated by means of the semiclassical expression 

<I -II- l|X_i(Af > =< M-iiMl) >= 



[/i3 sin "i? — /zi cos . 



(40) 



The components of the transition operator M.^, refer the 
the lab system. The expectation value is taken with the 
TAC configuration | >. In the second fine A^^ is ex- 
pressed by the components of the magnetic moment in 
the intrinsic frame. The reduced Ml-transition probabil- 
ity becomes 

< M-i{Ml) >2 . 



B{M1, I ^ I -I 
The spectroscopic magnetic moment is given by 

/ 



(41) 



=< II\^iz\II > = 

I 



1/2 



/ + 1/2 
[^1 sin + H3 cos ■&]. 



(42) 



The factor 



1+1/2 



is a quantal correction which is close 



to one for high spin. The components of the magnetic 
moment with respect to the PA are calculated by means 
of 

Ail = Miv(>^i,p + (?75-58 - 1) 5*1,,, - ?73.82 Si,n), 

= Miv( J3,p + ('75.58 - 1) S3,p - ?73.82 53,„), (43) 

where the components of the vectors of angular momen- 
tum J and of the spin S —< s > are the expectation 
values with the TAC configuration | >. The free spin 
magnetic moments are attenuated by a factor rj — 0.7. 
For mass other mass regions a somewhat different atten- 
uation may be taken, which needs not be the same for 
protons and neutrons. 

The intra band E2 - transition matrix elements are 
calculated by means of the semiclassical expression^ 

<I-2I- 2\M-2iE2)\II >=< M-2{E2) >= 
^ <Q'o> (sin i?)2 

+ i <Q'2 + 0-2 > (l + (cosz9)2)], (44) 
< / - 1/ - l\M-i{E2)\II >=< M-i{E2) > = 
sin i9cos??(w| < Qq > 



A , 



\<Q'2+ Q'-2 >) 



(45) 



For the treatment of the term X < N > see sect. 



IIJ 



"Refs. Q-g contain some unfortunate inconsistency be- 
tween the quadrupole moments of the quadrupole interac- 
tion and the electric transition matrix elements. These con- 
cern only the written formulae, the results of the calculations 
quoted are correct and consistent with the ones given here. 
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and the spectroscopic quadrupole moment by 



q^<ii\q[''^''>\ii> 



I 

/+2/3 



[BM] 



I 2eZ 
7+2/3 A 



<Q'o> ((cosi?)2 - i(sin 



< Q2 + 0-2 > (sin 1?)^ 



(46) 



We use the conventional definition of the static 



quadrupole moment as given in ref. |53|, which differs 
by a factor of 2 from our quadrupole moments in the lab 
frame. There is a similar quantal correction factor as for 
the magnetic moment. 

The reduced E2-transition probabilities are 



B{E2, I ^ I-2)=< M-2{E2) >^ 



and 



B{E2, / ^ / - 1) =< M-iiE2) >2 
The mixing ratio is 



6 = 



< M-i{E2) > 

< M-iiMl) >■ 



(47) 



(48) 



(49) 



The mass quadrupole moments consist two terms. The 
first one contains the microscopic expectation values 
< Qo>N and < Q'j_2 >Ar, where the subscript N indi- 
cates that only the AN — matrix elements of the 
quadrupole operator are taken. The second term takes 
care of the coupling between the oscillator shells. 

In the case of SCTAC the stretched coordinates are 
introduced to approximately take the coupling between 
the oscillator shells into account. The expectation values 
needed in eqs. (^, are the quadrupole moments 

in unstretched coordinates, which are given by 

<Q'o >= 

W2^_^^_^) 2 > 

D ^ CJ3 LJi UJ2 



I 1 / 4lJo 



't)<Q2 + Q'-2 >N 



(50) 



< 



,/K V, 



K Q-2 > = 
5(t + ^)<Q2 + 0'-2> 



N 



(51) 



where the semiclassical value < >n~ 1.2 Z A^^^^ fra 
is used. 

In the case of the PQTAC the coupling between the 
oscillator shells is neglected. This is a reasonable ap- 
proximation for the rotational response of the valence 
particles. However, when calculating electric quadrupole 
moments it cannot not be neglected, because it accounts 
for the polarization of the core by the valence nucleons. 



We describe the polarization by means of eqs ( pO| , pl| ), 
setting e = 5. This prescription satisfies the consistency 
condition that the deformations of the potential and the 
density should be the same . It corresponds to a po- 
larization charge close to 1, as estimated for the isoscalar 
quadrupole mode |]5^ . This choice of the polarization 
charge makes to PQTAC and the SCTAC as similar as 
possible. 

In the above described methods one could also use the 
proton part of the quadrupole moments instead of Z /A 
times the mass quadrupole moments. 



G. Quantization 

Due to leading quantal correction (cf. e. g. |^,^) one 
must associate the total angular momentum J calculated 
in TAC with I + 1/2, where / is the quantum number of 
the angular momentum. This prescription permits us 
to compare the TAC calculations with the experimental 
energies and the static moments. Genuine TAC solutions 
{d ^ 0°, 90°) represent A/ = 1 bands. In this case, the 
experimental rotational frequency uj is introduced by 

J ^ I, uj{I ^ I -I) = E{I) - E{I -I), (52) 

and the experimental Routhian by 

E'il ^ / - 1) = ^[E{I) + E{I - 1)] -u:{I^I- 1) J. 

(53) 

Here, the canonical relations ( |20| ) are approximated by 
quotients of finite differences. The data define a discrete 
sets of points J{oj) and E'(lu), which are connected by 
interpolation. If the axis of rotation coincides with one 
of the principal axes {d = 0°, 90°), states differing by 
two units of angular momentum arrange into a A/ = 2 
band of given signature a. In this case the frequency is 
calculated by 

J = I- 1/2, u{I^I-2)^ ![£;(/) - E{I - 2)], (54) 

and the experimental Routhian by 

E'{I I -2) = ]^[E{I) + E{I - 2)] -u){I ^ I - 2) J. 

(55) 

For the transition probabilities, J is associated with 
the mean value of /+ 1/2 of the transition, i. e. 

B{M1,I ^ I -1) = B{M1,J = I) (56) 
B{E2,I ^ I -1) = B{E2,J = I) (57) 
B{E2,I ^ I ~2) ^ B{E2,J = I -1/2), (58) 

where the rhs denotes the result of the TAC calculation 
taken for the indicated value of J. Another possibility 
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is to compare the experimental transition probabilities 
with the ones calculated at the experimental frequency 
of the transition (|5^,^). As long as the experimental 
and calculated functions J^lo) agree well, both ways will 
give about the same result. 

Of course, one may also use the relations (^J) and (|55| ) 
for a A/ — 1 band. Then, the two signature branches 
will lie nearly on top of each other if the discrete points 
are connected by smooth interpolation. This choice has 
the disadvantage that the distance between the discrete 
frequency points is doubled. It has the advantage to give 
smooth curves when the splitting between the two sig- 
nature branches gradually develops with increasing fre- 
quency. In such a case and ds^) should be used. 



H. Diabatic tracing 

The goal of the calculation is to describe a rotational 
band, which corresponds to the "same quasi particle con- 
figuration" for a set of increasing values of a;. This means 
that one should keep fixed the occupation of the quasi 
particle states with similar structure. Usually one band 
does not correspond to the same configuration if the quasi 
particle levels are labeled according to their energy, be- 
cause the quasi particle trajectories cross each other as 
functions of u and -d. In order to find the equilibrium an- 
gle, one has to calculate the functions Ji(i?) and J^ii})- 
This becomes very tedious if the configurations are as- 
signed manually by identifying the crossings from quasi 
particle diagrams like figs. || and ||. The task is greatly 
facilitated by tracing the structure of the quasi particle 
wave functions. The calculations are run changing or 
uj in finite steps. For a given grid point the overlaps of 
each quasi particle state with all states of the previous 
grid point are calculated. The pair with the maximal 
overlap continues one quasi particle level from the pre- 
vious to the present grid point. The pair with the next 
lower overlap continues the second quasi particle trajec- 
tory. This procedure is repeated until all quasi parti- 
cle trajectories are continued. For all the single particle 
and quasi particle diagrams shown in this paper the grid 
points are connected by means of this diabatic tracing. 

In a practical calculation, the configurations are as- 
signed manually for the first grid point in a loop. The fol- 
lowing strategy has turned out to be quite efficient: First 
a typical angle i^s is chosen and the quasi particle diagram 
e'{uji,'ds) is generated. The step size Aw = O.OSMeT^ 
has turned out to be a good choice. Configurations are 
assigned for a typical frequency. The occupation num- 
bers for the other grid points oji are found by means of 
the quasi particle diagrams or, if the crossing pattern is 
complex, using the tracing facility of the code. These 
occupation numbers are used to set the configurations in 
a i9-loop starting at uoi and i^s • The configurations of the 
other grid points in the loop are determined by means 
of diabatic tracing. Then, the code finds the orientation 



angle d for each by means of the self-consistency con- 
dition (UJ) and calculates the interesting quantities. The 
step size Ai9 = 5° has turned out to be a good choice. 
At which ds the loop is started depends on the type of 
the band and will be discussed below. 

Problems are encountered when the quasi particle lev- 
els do not cross sharply when •& oi lo are changing. If 
the grid point happens to be located in the middle of 
the region where the levels strongly mix and repel each 
other, the diabatic tracing does not always follow the 
desired structure. Such cases necessitate human inter- 
ference in order to continue the correct structure. One 
reruns the calculation with the complementary configu- 
ration and puts the parts with the correct configuration 
together. The grid point itself is problematic because the 
cranking model becomes a bad approximation due to the 
unphysical mixing of states with different angular mo- 
mentum. These problems have been investigated for the 
standard cranking model |5^. We restrict ourselves to 
the most simple solution advocated in We discard 
such grid points and bridge the crossing region by means 
of interpolation. 

The results of the diabatic tracing depend on the step 
size. It should not be too small. If the step size is much 
smaller than the mixing region, the procedure follows 
the levels adiabatically, i. e. it connects the levels (of 
the same parity) according to their energy. On the other 
hand it should not be too large in order to preserve a 
reasonable precision. As mentioned above, step sizes of 
M = 5° and Aw = 0.05 MeV have turned out to be 
good choices. 

For low K bands it is usually convenient to choose 
•ds — 85° for the manual assignment of the configurations. 
The reason is that with increasing lo the equilibrium an- 
gle do changes quickly from zero to values close to 90°. 
As seen in figs. || and ^, the number of avoided cross- 
ings, which cause problems, is small at low frequency. 
Therefore the diabatic tracing works well in most cases 
and permits calculating the interesting range of oj with- 
out human interference. The configuration assignment 
should not be done at 90°, where the signature is good 
and the levels are often degenerated . The configurations 
discussed in subsections HID III H are calculated by as- 
signing configurations at ds — 85°. 

For high K bands, 'do remains relatively small up to 
rather large values of u. Then, starting at i?s = 85° be- 
comes less efficient because the number of avoided cross- 
ings increases. A smaller value of -ds closer to the equilib- 
rium angle is preferable. For the configurations discussed 
in subsections [III we used "i?, — 45°. This choice has the 



disadvantage that one has to run the d loop two times, 
for -i? < i?s and d > ds- The choice — 0° has turned 
out to be quite efficient in other applications of TAG to 
high K bands. 

Diabatic tracing is also used when the other parame- 
ters of the mean field Hamiltonian are changed in order to 
solve the complete set of self-consistency equations. Ap- 
proaching the minimum on the multi dimensional surface 
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E'{uj, £, £4, 7, A) it is applied for each step in one of the 
parameters. 

I. Choice of the QQ coupling constant 

Using the PQTAC version, the couphng constant x of 
the QQ - interaction must be fixed. So far it has been 
adjusted such that the quadrupole deformations e and 
7 calculated for PAC solutions come as close as possi- 
ble to the ones obtained by means of the shell correction 
method, which has a considerable predictive potential 
concerning the nuclear shapes (cf. e.g. (s^]). The ad- 
justment has been carried out for selected nuclei. The 
QQ coupling constant scales with A~^l'^r'^^^, where rose 
is the oscillator length This scaling has been used 
to determine x neighboring nuclei. 

In the first TAG calculations Q the equilibrium shape 
was calculated for A — 170 using the standard shell cor- 
rection method at a; = 0. The calculation was repeated 
for PQTAC at the same deformation and = 0. The cou- 
pling constant x '^^s chosen such that the self-consistency 
equations (|3^) were fulfilled. This value of x was kept 
constant in the full TAG calculation for all values of u). 
The value x = 0.0174 MeV r~^^ was found. Scaling gives 
X = 91 MeV A~^l^r~^^ for the rare earth region. When 
SGTAG became available, it turned out that the results 
of PQTAG and SGTAG were nearly identical for the nu- 
clei around A = 170. 

Extrapolating by means of scaling gives 
X = 0.0133 MeV r"/^ for A = 200. This value (scaled 
locally) gives a good overall description of the magnetic 
dipole bands in the Pb- isotopes pT|-p^. A deformation 
of e « —0.11 is obtained. SGTAG gives larger deforma- 
tions of e « —0.15, which account less well for the data 
on the Pb - isotopes. 

Scaling of the rare earth value gives x — 
0.036 MeV r^t and 0.024 MeV r"* for A = 110 and 
140, respectively. A new adjustment of x. was car- 
ried out for ^^°Gd and ^^^Sm. The respective values 
X = 0.036 MeV r-'t and 0.022 MeV r"/^ were deter- 
mined by making equal the deformation obtained by 
means of PQTAG and the shell correction method for 
zero pairing at finite lo. The latter values (including local 
scaling) gave a good description of the magnetic dipole 
bands in a number of nuclides of the two regions p5| -p2[ . 
The data on electro-magnetic transitions in los.iOeTTOKg^ 
1^,^, which according to the calculations have a de- 
formation of e < 0.11, seem to point to a smaller de- 
formation than calculated. That is a smaller value of 
X leading to smaller deformation than predicted by the 
shell correction method appears to be more appropriate 
for i05406,i08gj-^ This is similar to the Pb-isotopes. 

The shell correction method accounts rather well for 
the overall tendencies of the shape, in particular for the 
well deformed nuclei. However, it is not obvious that for 
the small deformations encountered for magnetic rotation 



(typically |e| < 0.11) the shell correction method provides 
a reliable gauge for x- In such cases it seems preferable 
to fix the QQ coupling constant in a different way. Since 
X controls the quadrupole polarizability, one may adjust 
it to the static quadrupole moments of high spin states 
and the B{E2) values of transition between them, which 
are particular sensitive to the quadrupole polarizability. 
It seems promising to use this experimental information 
for a fine tuning of x- This approach, which is discussed 
in more detail in the review ||26| , is being investigated 
11- 

J. Approximate treatment of self-consistency 

The GHFB equations are a complex system of non- 
linear equations. TAG adds a new dimension to it, the 
orientation angle -d. The fully self-consistent solution of 
the equations becomes rather tedious, in particular if one 
tries to describe several non yrast bands. The success of 
the GSM shows that for a first analysis of the ex- 
citation spectrum it is often sufficient, even preferable, 
to keep fixed the parameters of the mean field. For se- 
lected bands, they may be determined self-consistently 
in subsequent calculations if one is interested in specific 
properties. But very often the additional effort docs not 
pay off the gain in insight. We shall follow the GSM ap- 
proach and carry out the calculations assuming that the 
deformation and the parameters of the pair field, A and 
A do not depend on the rotational frequency uj. Only the 
tilt angle, which usually strongly changes, is determined 
by means of (|^ for each value of w. 

For the well deformed nuclei considered in this paper 
the deformation changes turn out to be moderate. They 
are negligible for the more qualitative comparison with 
the data which we are aiming at. 

The approximation of a constant pair field needs a 
more careful discussion. The original assumption of the 
GSM to keep A at 80% of the experimental odd - even 
mass difference Aoe becomes problematic, because mod- 
ern data reach rotational frequencies where the static pair 
field disappears | |59[ | . Since the transition to the unpaired 
state may substantially change rotational response, self- 
consistency must be taken into account at least in some 
rough way. We found the following compromise between 
accuracy and effort quite satisfying. The TAG calcula- 
tions are carried out at a few values of A, which do not 
depend on uj. The total Routhians E'{ui, A) are plot- 
ted. At each frequency one can easily choose the best 
A- value as the one that has the lowest value of E' . The 
upper panel of fig. ^ shows as an example the yrast 
band of "-^Hf. The values A„ = 0.95, 0.69, MeV 
and Ap = 1.05, 0.75, MeV are used. The first point 
corresponds to the self-consistent ground state value 
and the corresponding curve is lowest for small uj. At 
u! — 0.24 MeV the curve with the reduced values of 
A„ and Ap takes over. Within the considered frequency 
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range the unpaired solution cannot compete, though the 
neutron correlation energy is rather small. 

For paired configurations the proper Routhian is 
E'{uj,X) + XN, where N is the exact particle number. 
The term XN compensates the Lagrangian multiplier in- 
troduced in (H). However, exact compensation appears 
only if the self-consistency condition (||) for A is fulfilled. 
If A is approximately determined one can proceed in a 
similar way as for A by plotting E'{uj, A) -I- AA^. Since 

— i?'(c^,A) = -<iV>, — <iV> >0, (59) 

the Routhian E'{lo, X) -I- AiV has a maximum at the self- 
consistent value of A. Accordingly, TAG calculations are 
carried out for a few values of A, which do not depend on 
Lo and E'{lu, A) -I- AA^ is plotted. The highest curve corre- 
sponds to the best value of A. The lower panel of fig. ^ 
shows the three points A„ = 49.0, 49.15, 49.3 MeV. The 
arrows indicate the frequencies where the self-consistency 
condition (§) is fulfilled. For A = 49.0 MeV, the devia- 
tion in particle number is about 2 at w = 0.5 MeV. The 
upper envelop of the curves represents the best choice of 
A and A within the restricted set of grid points investi- 
gated. For uj > 0.3 MeV, it behaves very similar to the 
unpaired curve. The small correlation energy of 0.1 - 0.2 
MeV indicates weak static pairing. We will show this 
optimized Routhian of the yrast sequence in the figures 
as a reference. As long as the values of A are the same for 
all configurations one may leave away the term AA'^. It 
is however needed to correctly calculate the relative po- 
sition of configurations with different A or of paired and 
unpaired configurations . 

This method is quite useful because it is simple and 
it can easily be made as accurate as needed by adding 
more A and A values. At each stage one has a clear idea 
of the remaining error of the energy. The simplest vari- 
ant of considering only A « O.SAoe and and choosing 
Xp and A„ such that the particle numbers are right for 
a; = turns out to be sufficient for a first orientation. It 
shows the pair correlation energy directly. The discus- 
sion of pairing will be restricted to this minimal variant. 
All figures showing total Routhians display the quantity 
E' + XpZ + XnN. In order to keep the figures simple, 
the ordinate is labeled with E' only. The energy of the 
ground state of the nucleus, which is not the concern of 
this paper, is not calculated correctly. In all figures, only 
the Routhians relative to the ground state energy are of 
relevance. 



K. Reading the quasi particle diagrams 

Figs. H - H show the single particle Routhians e^{uj, 
as functions of the frequency u and of the tilt angle ?9. 
In order to demonstrate change of the particle response 
to with the magnitude and orientation of uj two different 
frequencies are presented for each kind of particles. The 



figures with intermediate frequency are relevant for the 
present day high spin data. The high frequency figures 
show territory yet to be explored. The side panels of each 
figure are added for helping the reader to connect to new 
middle panel with the familiar single particle Routhians 
for rotation about the PA axes. 

The slope of the trajectories gives negative projection 
of the quasi particle angular momentum on the u - axis, 

de' 

-Q^ = = -(ji,jsint?-h ja.jcost?), (60) 

and its perpendicular component, 
de' 

= -^j± = cos?? - js^, sin-d), (61) 

where ji^i and js^i are the expectation values of the an- 
gular momentum in the single-particle or quasi particle 
state i. 

For ?? = 0°, the cranking term L0J3 commutes with the 
axial symmetric deformed potential and the projection 
of the angular momentum on the symmetry axis K is a 
good quantum number. In this case, the states coincide 
with the non-rotating Nilsson states, which are indicated 
by the labels in the figure. For •& = 90°, the signature 
a is a good quantum number which is also indicated in 
the figure. As the signature operator e^"'^^ and Js do 
not commute, there is a transition from one to the other 
type of symmetry when i? changes from 0° to 90°. 

Discussing the features of the quasi particle diagrams, 
we will refer to the three types of coupling schemes that 
appear as a consequence of the competition between the 
deformed potential, the inertial forces and the pair cor- 
relations. They are discussed in | |60[ |. Let us start with 
moderate frequencies w < 0.3 MeV , which are illustrated 
in figs. H, and 0. 

The normal parity states with X > 5/2 obey the defor- 
mation aligned coupling (DAL) scheme. These orbitals 
are strongly coupled to the deformed potential. In the 
e'(i?) plot, they are recognized as the pairs of trajecto- 
ries, which branch at d = 90°. They have a small com- 
ponent ji^i but a large component j^^i « K. They ap- 
proximately behave like —Kcos'd in an extended region. 
Near d = 90° there is the very narrow transition region 
from good signature to almost good K, where the slope 
changes from zero to approximately Ksiwd. The region 
is too narrow to be discerned in the figure, where it looks 
like a kink. 

There are pairs of parallel trajectories in the e'(i?) plot, 
originating from the [521] 1/2 neutron and [411] 1/2 pro- 
ton orbitals. These are pseudo spin singlets. A discus- 
sion of the pseudo spin symmetry in deformed poten- 
tials is given in [Q. The projection A3 of pseudo orbital 
momentum is zero. Thus, the pseudo spin is decoupled 
from the orbital motion. It only reacts to the crank- 
ing term — cJ • S, where S is the pseudo spin. The two 
parallel, nearly horizontal trajectories with the distance 
Acpso = w correspond to the pseudo spin being aligned 
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or anti aligned with the rotational axis uj. The pseudo 
spin vector follows the tilt of w, remaining parallel to it. 
Since the pseudo orbital momentum remains small, the 
two trajectories are almost horizontal. The signature a 
is gradually lost when the pseudo spin vector tilts away 
from the 1 - axis. 

The states with the highest K values of the /iii/2 and 
*i3/2 intruder orbitals obey the DAL coupling. The states 
with lower K have an extended region around 90° where 
the Routhians are relatively flat functions of i), what 
means that j± is small. In this region the orbitals are 
rotational aligned (RAL), precessing around the rota- 
tional axis uj, where the precession cone follows the tilt 
of the axis. The signature is gradually lost when uj tilts 
away from the 1 - axis. With decreasing d, they make a 
quasi crossing with other members of the same intruder 
orbital. These crossings mark the transition to the de- 
formation aligned (DAL) coupling, which is shows up as 
the —Kcosd behavior. 

Fig. ^ shows the quasi neutron energies, which are rele- 
vant when the pair correlations are important. What has 
been said about the single particle Routhians also applies 
to the quasi particle Routhians. As a new type, the Fermi 
aligned (FAL) couphng Q appears. It is reahzed by the 
lowest ii3/2 trajectory, denoted by A. The FAL couphng 
appears at some distance from 90°. It corresponds to a 
substantial component as well as to a substantial ji. 
It is most favored at the minimum of e^(??) a,t — 38°, 
where j± = 0, that is j\\uj. With i9 — > the non rotating 
quasi particle state is approached, i. e. — > K and 
ji 0, corresponding to the maximum. Overall, the 
lowest ii3/2 trajectory A is rather flat, indicating that 
the orientation of j does never too strongly deviate from 
Lu. At larger uj, where the negative and positive quasi 
particle states strongly interact with each other, a com- 
plex pattern of avoided crossings emerges, which we have 
not found a simple interpretation for. 

The high frequency regime is illustrated in figs. ^ and 
^. It is characterized by many avoided crossing between 
the orbitals. This indicates the progressive dissolution of 
the approximate conservation of the K quantum number 
for the DAL orbitals. 



L. Relation to the PAC treatment of high-7f bands 

Bands with a finite value of K have been studied by 
means of the standard PAC scheme using the following 
prescription p7|| , which may be considered as an approxi- 
mation to TAG. A fixed K value is ascribed to each band, 
which is the spin value at the band head. It is taken 
from experiment or calculated by means of the cranking 
model choosing lu parallel to the symmetry axis. It is as- 
sumed that J3 = K, independent of uj and Ji =< Ji >. 
The configuration | > is generated by the quasi parti- 
cle Routhian (p2h assuming — 90°. Only the reduced 



cranking term — wi Ji appears, where uji is the 1 - com- 
ponent of the angular velocity. 

With these assumptions TAG goes over into the GSM 
1^ scheme: The constraint ( |l7|) becomes 



J = 1+1/2= Jjl + J| = V< Ji +K'^: (62) 



Solving for < Ji >, the standard cranking constraint 



< Jl >= Jlv/(/+l/2)2-if2 



(63) 



to fix uJi is obtained. In the GSM one uses uji as the inde- 
pendent variable. Experimental values of uji are derived 
by means of the expression 



E{I)-E{I-2) 
Ji(/)- Ji(/-2)' 



(64) 



with Jl being the rhs. of expression (pSj). The TAG 
condition J 1 1 implies 



UJl 



Jl _ - 1/2)2 - if 2 



J 



2(7-1/2) 



{E{I)-E{I-2)), (65) 



where the expression (|5J) for the experimental frequency 
is used. The uji values obtained by expressions ( |64|) and 
( |65| ) almost coincide, except near the band head. As 
demonstrated in the model study expression ( |65| ) 
reproduces the quantal results slighly better. 

With the fixed- ivT assumption the expressions for the 
electro-magnetic matrix elements given in sect. II F be- 
come the ones of the semiclassical vector model of ref. 
psf . For the magnetic moments, the vector model ad- 
ditionally assumes that each quasi particle has a fixed 
value of j3_i, which is uj independent and given by the 
Ki of the Nilsson label. The individual ji^i values of the 
quasi particles are either calculated or extracted from 
differences between the experimental curves Ji(wi) ( the 
quasi particle alignments ). The magnetic moments are 
approximated by 



(66) 



where the gyromagnetic ratio is either calculated by 
means of the Nilsson model or taken from experi- 
ment. In TAG expression ( p| , [12D the components of the 
magnetic moments are calculated, attenuating the free 
spin magnetic moment of the proton and neutron by a 
constant factor. 

In ref. an additional term is introduced which per- 
mits the calculation of the signature dependence of the 
B{M\) values for the case that the J3 component is gen- 
erated by only one quasi particle . It is not expected 
that this correction is also applicable if J3 is generated 
by many quasi particles , whereas the vector model ]6^ ] 
without the signature term also applies to this more gen- 
eral case. 

As discussed, the standard GSM becomes a good ap- 
proximation of TAG if for the active quasi particles 
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i) j3^i{u!, 'd) can be approximated by the constant Ki and 

ii) ji.i{uj, ^) can be approximated by 90°). 
Then 



s^L^,^) «eKc^i,90°)-c^3^, 



(67) 



where ui = wsini? and 1^3 = uj cost} (cf. eqs. ( [3^ ) and 
(^)). The assumption J3 w X is justified and the tilt 
angle is given by cos d = K/ Jj, which leads to the ex- 
pressions of the vector model . 

Fig. ^ shows the angular momentum components 
ji{ijj,-d) and j^{uj,'d) of some representative quasi par- 
ticles. They are compared with the CSM values 
ji,i(wi,90°) and Ki, respectively. The DAL quasi neu- 
tron E obeys i) and ii) rather well. The orbital G shows 
the behavior ji « i sini^ and js ~ i cosz^, which is char- 
acteristic for the pseudo spin singlet. This is at variance 
with i) and ii). Although the contribution to the an- 
gular momentum is small, the characteristic spacing of 



uj between the two pseudo spin partners is not 



obtained when -d is substantially below 90°. The intruder 
orbital A shows in the range 40° < ■& < 60° the typical 
FAL behavior, which is fairly well reproduced by the ap- 
proximations i) and ii) . For larger values of § it changes 
gradually into a state with a good signature. The transi- 
tion is accompanied by changes of ji and ja which are at 
variance with i) and ii) . For smaller values of d orbital A 
keeps its FAL character and remains close to the approx- 
imations i) and ii). The orbital B changes dramatically 
when 1? decreases from 90°, because there is the quasi 
crossing with the down sloping orbital C (cf. fig. ^ . For 
values of m smaller than displayed this crossing is rather 
sharp. One can follow the FAL branch of B, which is 
well approximated by i) and ii), below the crossing. For 
larger values of oj the two orbitals strongly mix and a new 
•d dependence emerges, which is shown fig. |ll|. Obvi- 
ously such changes of the quasi particle structure cannot 
be described by means of the traditional CSM treatment 
basing on the assumptions i) and ii). 

The use of wi as the rotational parameter in the CSM 
has the advantage that all configurations can be con- 
structed from one and the same quasi particle diagram 
ei{u!i, 90°), which is a great simplification and has lead to 
the popularity of this approach. Another pleasant feature 
is that the signature splitting appears in a gradual way. 
The disadvantage is that it is only an approximation to 
the TAC mean field solution, the latter being completely 
self-consistent and more accurate when d substantially 
deviates from 90° ||6^ . Sometimes the differences are only 
of quantitative nature, but there are many cases where 
qualitatively different results are obtained. Magnetic Ro- 
tation of weakly deformed nuclei |]lT|-p^ is a conspicuous 
example, which will not be discussed in this paper. In 
the following discussion of examples we shall point out 
the differences between TAC the standard CSM. 

High- if bands, which are in experiment near yrast, ap- 
pear in the CSM as relatively high lying configurations, 
embedded into the back ground of many configurations 



with low K. In TAC they are low lying configurations. 
The reason can be seen in (|^. CSM uses e-(wi,90°), 
which is shifted up by to^Ki with respect to e[{uj, d) used 
in TAC. It is also noted that only the TAC mean field so- 
lutions can be improved by means of RPA corrections in 
a systematic way. The PAC configurations corresponding 
to a finite K are instable. 

It is quite common to present the experimental branch- 



ing ratios 



as effective values of the ratio 



which would determine the branching ratio if the strong 
coupling limit was valid |5^. This popular way of rep- 
resenting the data has the advantage that the ratio be- 
comes constant when approaching the strong coupling 
limit. One may convert the ratios of "^''^'^^j* calculated 



by means of TAC into effective ratios 
tinent relation 



Qo 



9K - 9R 



Qo 



5(/ 



^ - K'^)B{Ml) y 
\Q,PB{E2) ) 



The per- 



(68) 



is obtained from the expressions in section [I F by mak- 
ing the assumption of strong coupling, J3 = if, /i3 = 
{9k — 9r)K and /ii = 0. The square root of the branch- 

and the inverse 



ing ratio becomes the product of 

of the geometric factor on the right hand side of (|6S^ 
In order to avoid any misunderstanding it is noted that 
( |68| ) is just a way to present the results of the exact TAC 
calculations, which do not make any strong coupling ap- 
proximation. 



III. MULTI-QUASI PARTICLE 
CONFIGURATIONS NEAR "*HF 

This section will explain how to construct multi quasi 
particle configurations in the TAC scheme and how to in- 
terpret them as rotational bands. The nuclides ^^^'^^'''Hf, 
^^^Ta and ^^^Lu serve as examples. The SCTAC scheme 
is used for the calculations. In order to simplify the 
discussion, the same deformation parameters e — 0.258, 
£4 = 0.034 and 7 = 0° are assumed for all the nuclides 
considered. The set represents an average of equilibrium 
shapes calculated for several configurations and frequen- 
cies Lo. The actual values scatter within the interval 
0.25 < e < 0.27, but the differences in deformation do not 
change the discussed quantities in a substantial way. We 
consider both the cases of no pairing and a constant pair 
field. For the case of finite pairing we use the prescription 
of the CSM 1^^, which has turned out to give very rea- 
sonable description of multi-quasi particle bands in the 
traditional PAC scheme. Accordingly, A„ = 0.69 MeV 
and Ap = 0.75 MeV, which is 80% of the experimen- 
tal even -odd mass difference. The chemical potentials 
A„ and Ap are fixed to the values that give the correct 
particle numbers for the ground state (configuration [0] 
at = 0). This scenario provides a good description 
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of configurations up to two excited quasi particles and a 
frequency o f abou t 0.35 MeV. It will be discussed in the 
subsections III D - [II H . For higher frequency and more 



excited quasi particles we take zero neutron pairing into 
consideration. Self-consist ency for A is invoked along 
the lines described in sect. [IJ in order to decide where 



the transition to zero pairing is located. Subsection HI J 
discusses this regime. 



A. Construction of multi quasi particle 
configurations 

For zero pairing the configurations are generated by 
filling up the lowest Z and N single particle levels and 
then making particle-hole excitations. 

In the case of pairing the configurations are con- 
structed from the quasi particle Routhians. The quasi 
particle spectrum is symmetric with respect to e' = and 
the double dimensional occupation scheme, discussed for 
the PAC solutions in ref |4^, is applied: If a quasi parti- 
cle state i is occupied, its conjugate partner i+ must be 
free. In contrast to the PAC case, the conjugate states 
in general do not have opposite signature, which is only 
for § = 90° a good quantum number. 

Diabatic tracing turns out very practical for identifying 
the conjugate states. They always cross sharply because 
they are orthogonal. As in the PAC scheme, one must be 
careful in choosing the right particle number parity when 
quasi particle trajectories cross the zero line. The most 
simple way is to start at sufficiently low to, where there 
is still a gap between the positive and negative solutions. 
There it is clear how to excite an odd or even number 
of quasi particles . Keeping the occupation by diabatic 
tracing, the particle number parity of the configuration 
is conserved. 

In order to efficiently label the configurations a com- 
pact notation which indicates the quasi particle composi- 
tion is desirable. We follow the well-tried practice of the 
CSM assigning letters to the quasi particle trajectories 
and quoting the excited quasi particles in parenthesis. 
The letters A, B, C, D denote positive parity quasi neu- 
trons , E, F, G, H, ... negative parity quasi neutrons , 
a, b, c, d positive parity quasi protons and e, f, g, h 
negative parity quasi protons . 

The letter code becomes to some extend ambiguous 
when the structure of the quasi particles strongly changes 
with the frequency lu and the tilt angle i?. The positive 
parity ii3/2 orbitals are most susceptible to the inertial 
forces. Fig. || shows the complex pattern of quasi particle 
trajectories, which strongly interact with each other and 
interchange their character as functions of uj and z?. An 
example are the orbitals B and C in fig. |3. As discussed 
already in sect. 



II L , they quasi-cross each other near 



1? = 60°. For LU — 0.2 MeV, (not shown) the crossing is 
still rather sharp and it would be natural to follow each 
trajectory diabatically, i. e. for i? < 60° to call the upper 



trajectory B and the lower one C. For ui = 0.4 MeV (see 
fig. pj]) they interchange their character very gradually. 
Now it is more natural to call the lower trajectory B 
and the higher C throughout the mixing region. Fig. 
[To| shows the quasi neutron trajectories as functions of 
uj. For d = 45° (lower panel) the crossings are rather 
sharp for most trajectories. Here it is natural to keep the 
labels in a diabatic way, as indicated. For ■& — 90° the 
crossings between the «i3/2 trajectories are much softer 
and the question arises of how to label them after the 
first quasi crossing. The suggested labefing tries to follow 
the quasi particle trajectories both in the uj and the 
direction such that the structural change is as gradual 
as possible. It connects the two diagrams e-(w, d = 45°) 
and e^{uj, z9 — 90°) the most natural way. via the i) degree 
of freedom at high frequency (cf. figs. |^ and |ri|). This 
implies that the smooth crossing between A and B+ as 
function oi uj at d — 90° must be treated diabatically. 

It should be pointed out that the suggested labeling 
is a compromise. As discussed above, for uj = 0.2 MeV 
the quasi neutrons B and C cross sharply as functions 
of d. In the adopted labeling the lower of the two levels 
is B and the higher C. It would be more natural to fol- 
low the structures in 'd direction diabatically through the 
crossing. But a relabeling that accounts for this leads to 
problems a high uj, where the adopted labeling is most 
natural. The difficulty to label the strongly interacting 
quasi particle trajectories in a simple way has a topolog- 
ical origin, which can be best understood if one follows 
in one of the quasi particle diagram H, ^, ||, ^ a trajec- 
tory on the (w,i9)-path: (0,90°) ^ (0.3 MeV,90°) -> 
(0.3 MeV, 0°) (0, 0°). For weakly interacting trajecto- 
ries, as most with normal parity, one returns to the same 
quasi particle . For the intruder trajectories, as ii3/2 and 
hii/2, this is not always the case. 

Trying to keep the notation as simple as possible we 
shall assign the low-w composition to a configuration and 
shall not change it when a crossing is encountered. The 
structural change can be figured out from the quasi par- 
ticle diagrams. 



B. Elimination of spurious states 

Each configuration with an equilibrium angle -do 7^ 90° 
is associated with a A/ = 1 rotational band (TAC so- 
lution), whereas each configuration with an equilibrium 
angle lio = 90° is associated with a A/ = 2 rotational 
band (PAC solution). Of course, the number of quantal 
states cannot abruptly double when the equilibrium an- 
gle moves away from 90°. Hence, one has to be careful 
in avoiding spurious states. This problem was studied 
in ref. | |6^ for the model system of one and two quasi 
particles coupled to a rotor. An elimination scheme has 
been suggested which is based on the following principle: 
The number of TAC configurations must be the same as 
the number of PAC configurations at d — 90°, where 
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they emanate from. This means, for each TAG mini- 
mum, which is interpreted as a AI = 1 band composed 
of two AI — 2 sequences, one has to discard one con- 
figuration. One finds this spurious state most easily by 
taking into account that the function E'(i!)) is symmetric 
with respect to 90°. If one configuration (J3 Ki K > ) 
has a minimum at do its mirror image [K < ) has the 
minimum at 90° — do- Tracing the function E'{d) of the 
mirror image diabatically through d = 90°, one arrives 
at the spurious configuration that must be discarded. 

The tilt angle do increases with lo. When it approaches 
90° one has to switch from the TAG to the PAG inter- 
pretation. This results in a discontinuity of the function 
E'{uj) for the unfavored (upper) signature branch: In the 
TAG scheme it is is degenerate with the favored branch 
whereas in the PAG scheme it is the discarded configu- 
ration , which is now taken into account. The angle for 
switching from TAG to PAG is to some extend arbitrary. 
We have found it reasonable to use the PAG interpre- 
tation when the equilibrium angle do > 80°. If several 
quasi particles combine into high-/C and low- if configu- 
rations it is important to switch from TAG to PAG for 
both the high- and the low- if configurations at the same 
ui. As discussed in |^ and in sect. [II G for a concrete 
example, changing to PAG only for a part of the config- 
urations results in highly nonorthogonal states. 

As an example, let us consider the most simple case 
of the one quasi neutron configurations denoted by [E] 
and [F] in fig. which represent, respectively, the two 
branches « 5/2 and -5/2 of the DAL orbital [512]5/2. 
For ui = 0.2 and 0.3 MeV, [E] has a minimum below 80", 
which is interpreted as a AI = 1 band. The diabatic 
continuation of [F] becomes the mirror image of [E] for 
d > 90°. Hence, [F] is spurious and must be discarded. 
The kink - like minimum of the upper branch of [5 12] 5/2 
must also be disregarded. 

The elimination rules are somewhat differently formu- 
lated in ref . . The reader might find this complemen- 
tary formulation instructive. The proposed scheme has 
been tested for the model system of one and two quasi 
particles coupled to a rotor . No spurious states have 
been found in the low lying spectrum after applying the 
elimination rules. 



C. Band heads 

Generally, a band is a quasi particle configuration 
whose angular momentum increases with the rotational 
frequency uj. Its structure changes gradually with w, such 
that it remains similar for adjacent quantal states of the 
band. This is a natural definition which permits calculat- 
ing both the start and the termination of a band. Since 
we restrict ourselves to well deformed nuclei, we shall 
discuss only the start in this paper. 

Fig. |l^ illustrates how the configuration [E] in ^^^Hf 
starts. The function E'{uj,d) has a minimum at = 



(and 180°) for uj below the band head frequency of 
ujfi — 0.08 MeV. In this range of u> the band has not yet 
started, because angular momentum does not depend on 
Lu, being J = K. The band actually starts at loh when 
the equilibrium value do becomes finite, i. e. when the 
minimum of E'(d) at d — 0° turns over into a maximum 
and there appears a minimum at do > 0°. That is, the 
frequency ujh where 



d^E'{LUh,d) 

dd^ 



i3=0 = 







(69) 



has the physical meaning of the rotational frequency of 
the band head. Fig. |lj shows the equilibrium angles do 
for several one quasi neutron bands. The bands heads lie 
where do bifurcates from the zero line. 

The experimental band head frequency is the energy 
of the first transition I = K+1^I = K. It should be 
compared with the frequency where TAG gives J(w) = 
K + 1, which is somewhat larger than w/j. In some of 
the figures (e. g. |l^) this frequency is indicated by a fat 
dot. In most of the figures the calculated curves start 
with the first grid point for which do > 0, i. e. Uh is only 
determined with the accuracy of Auj — 0.05 MeV, which 
is the step used in the calculations. 

One may distinguish between strong coupling behavior 
and more complex response near the band head. Strong 
coupling behavior corresponds to J3 = if = const, and 
Ji = LiJiJ, where J is the moment of inertia of the col- 
lective rotation. In this case one has 

K K 
J = ujJ, = arccos( — ) for uj > uJh ^ —. (70) 
i-^J J 

Axial nuclei are close to the strong coupling limit near the 
band head if only DAL quasi particles are excited. The 
configuration [E] illustrated in figs. |l^ - ^is of his type. 
At a first glance, one might expect that the TAG approx- 
imation, which treats the orientation angle in a static 
way, becomes a bad approximation near the band head, 
because the E'{d) is very flat there. The model studies 
p2[ demonstrated that this is not the case. In fact, the 
wave function becomes narrow at the band head, because 
J3 approaches the good quantum number K. One may 
interpret this as follows. The mass parameter associated 
with the zero point motion in d increases faster than the 
curvature near the minimum, d'^E'(uih,d)/dd'^ |i?=iJo- 

A more complex situation is encountered when one or 
more quasi particles easily align with the rotational axis. 
Configuration [A] in fig. |l^ is an example. The band 
starts significantly earlier than expected for a strongly 
coupled 7/2 band with a jump of do. For cases like this 
ref. [ |62[ found that TAG approximates the quantal par- 
ticle rotor calculation less well near the band head, but 
becomes again a very good approximation for higher fre- 
quency. 
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D. The zero quasi particle configuration 



E. One quasi neutron configurations 



In ^^''Hf, the lowest configuration at i? = 90° is the 
vacuum [0] with all negative levels occupied, which is 
the ground (g-) configuration at low lo. Around lo = 
0.30 AleV, the neutron system gradually changes into 
s-configuration, which is seen in figs. ^ and (upper 
panel) as the quasi crossings of trajectories A with B+ 
and B with A+ (AB crossing). Near these crossings the 
1? dependence of the trajectories is complicated. We shal l 
return to the interpretation of this region in sect. IIIH. 



First, let us discuss the proton system, which does not 
have such a crossing in the considered frequency interval. 
At ■& = 90°, the configuration [0] has the character of 
the g-band. It keeps this character when d decreases, 
provided the occupation is followed diabatically, i. e. the 
crossings between the tt = + trajectories at z9 = 22° 
and the tt — — trajectories at "ff — 8° are ignored. It 
becomes the ground state for z9 = 0, because the wave 
function does not depend on uj for this orientation. The 
ground state is not the lowest configuration at i9 = 0, 
because a number of quasi particles have crossed the zero 
line and crossed each other. This example demonstrates 
the advantage of diabatic tracing, which automatically 
finds [0] when starting from either small uj and = 0° 
or 1? close to 90°, where [0] is the lowest configuration 
. The neutron system has an analogous structure for 
UJ < 0.25 MeV, where where [0] has the character of 
the g-band. Fig. ^ shows the quasi neutron trajectories 
at w = 0.3 MeV. It is seen that the configuration [0], 
which has a mixed g- and s- character d — 90°, becomes 
the ground state at i9 = 0°, where it is no longer the 
lowest configuration . 

Fig. ^ shows the total Routhian E'{uj = 0.2 MeV,d) 
of the combined proton and neutron configurations [0]. 
Its 1? dependence refiects the g-band character: The an- 
gular momentum is collective, i. e. 



{lu sini^)^ 



J. 



(71) 



The minimum lies a 1} — 90°, the signature is a = 0, cor- 
responding to the even spin g -band. For lo = 0.3 MeV 
the level repulsion near between A and B+ modifies the 
slightly the i? dependence of the total Routhian. 

In order to make a first qualitative estimate of the tilt 
angle for multi quasi particle configurations one has to 
add this zero quasi particle Routhian to the sum of the 
Routhians e^(?9) of the excited quasi particles , which can 
be taken from the quasi particle diagrams. 

For u! > 0.35 MeV the quasi particle vacuum has the 
character of the s-configuration a,t -d — 90°. With de- 
creasing it changes into the t-configuration which be- 
comes the K"" = 8+ configuration [7/2+,9/2l a.t ■& = 0°. 
We s hall di scuss these changes and their consequences in 
sect. IIIH, together with the two quasi neutron excita- 
tions of positive parity. 



They are generated by adding one quasi neutron to the 
configuration [0]. Fig. |l^ shows their total Routhians 

The configurations [G] and [H] have -do = 90° for all uj. 
They are interpreted as A J = 2 bands. They represent 
the two signatures a = ±1/2 of the pseudo spin singlet 
[521] 1/2. Since the pseudo spin is decoupled (cf. sect. 
[TkI ), E'{d) ^ E[^^{d)+const±uj/2 and = 90°. 

For u! = 0.2 and 0.3 MeV the configurations [A] and 
[E] have minima at 90° > iJo > 0°. They are interpreted 
as A/ = 1 bands {K'' = 7/2+ and 5/2"). The con- 
figurations [B] and [F] are the continuations of [A] and 
[E] reflected through i? = 90°. Accordingly they are dis- 
carded as spurious states together with the kink at 90°. 
Then the condition is satisfied, that the number of states 
is the same as for the PAC interpretation at = 90°. 
For LU = 0.4 MeV the minima of [A] and [E] have moved 
above 80°. Now we change to the PAC interpretation 
and refer to the calculations at •& = 90°. Both [A] and 
[B] are interpreted as A/ = 2 bands. They form the sig- 
nature pair (tt, a) = (+, ± 1/2). The configurations [E] 
and [F] represent the two A/ = 2 bands combining to 
the signature pair (— , ± 1/2). 

Fig. 16 shows the calculated total Routhians. The 



change from the TAG interpretation to the PAG is seen as 
the sudden onset of the signature splitting. The width of 
the transition region is determined by the size of the cal- 
culation grid in lu, which is 0.05 MeV and does not bear 
any physical relevance. As discussed above, the TAG ap- 
proach is not able to describe the smooth transition from 
broken to conserved signature symmetry. In order to de- 
scribe the gradual onset of the signature splitting one has 
to go beyond the pure mean field theory ||7|, ^ , p6|| . 

Since the orbital E obeys the DAL coupling, the con- 
figuration [E] is expected to be close to the strong cou- 
pling limit. Fig. ^ shows that the 5/2" band starts 
at LUk. — 0.09 AleV near the strong couphng estimate 
Kj J . Also for higher uj the tilt angle i?o remains close 
to the strong coupling value. The 7/2+ band starts at 
tj^j — 0.02 MeV , below the 5/2~ band [E] and much be- 
low the strong coupling estimate for the K = 7 /2 band. 
This indicates a substantial deviation from strong cou- 
pling. In fact, fig. |l^ shows that -do jumps to a finite 
value at a low frequency, corresponding to the rapid tran- 
sition from the DAL to FAL coupling at the band head. 

Fig. |l6| also shows the experimental Routhians [|66| . 
The relative position of the Routhians as well as their 
slopes (i. e. the angular momentum ) are reasonably 
well reproduced by the calculation. The frequency of 
the first transition is well described too. In particular 
the low value of LUh for the 7/2+ band indicates that the 
FAL coupling is seen in the experiment. 

In the TAG calculation the configuration [G] starts at 
LUh = 0.19 MeV with J = 6.1 and = 60° (cf. fig. 
It represents the K = 9/2+ orbital. The configuration 
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[D] is discarded because it is the continuation of the mir- 
ror image of [C]. The minimum rapidly moves towards 
90°, due to the strong admixture of 113/2 components 
with low K. After uj = 0.25 MeV one has to change 
to the PAC interpretation, where both [C] and [D] are 
interpreted as the signature pair (+, ± 1/2). Experi- 
mentally, only the transition / — 9/2 — > 11/2 is seen at 
UJ = 0.153 MeV, which is lower than the calculated value 

of LVh- 

The panel uj = 0.2 MeV in fig. |l| shows that the 
minimum of [C] is very shallow. For oj = 0.19 MeV it 
becomes a shoulder. In contrast to the experiment, there 
is no ?? > solution for lower frequency. Here, a limita- 
tion of the TAG approximation is encountered. The tilt 
angle i^o is found in a static way by searching for the min- 
imum of the Routhian This is an approximation 
to studying the dynamics of the "d degree of freedom. The 
static TAG treatment is expected to give good results as 
long as there exists a certain convex region around t9o- 
Then 'd will execute symmetric oscillations and averaging 
over them will result in values close to the ones for t?o. 
The model studies in ref. have demonstrated this for 
the lowest configurations. It is clear for a curve like [G] 
that averaging over the the collective wave function in 
needs not to give values close to the ones obtained for t9o, 
in particular, when the minimum has become a shoulder. 
Then the dynamics of d must be explicitly calculated. 
Refs. [[zt ^Sp^ have addressed this problem in the frame 
work of the Generator Coordinate Method. 

Situations like the discussed one become more likely 
if one considers excited configurations . As seen in fig. 
[1^ , the flat behavior of [C] may be thought as the con- 
sequence of interaction (repulsion) with the configura- 
tions below and above. This problem is not special to 
the orientation degree of freedom. Analog restrictions of 
the static HFB approximation are encountered when it 
is used to calculate the shape of excited configurations. 

Since the J3 < 4.5 for all the one quasi neutron con- 
figurations the tilt angle 'do rapidly increases with the 
frequency. As seen in fig. |lj, t9o > 60° for u > 0.2 MeV. 
Accordingly, the quasi particle trajectories become sim- 
ilar to the ones aX 1) = 90°. One recognizes the familiar 
GSM pattern of band crossings. The tt = — bands show 
the AB crossing and the tt = + bands the delayed BG 
crossing, because AB is blocked (cf. e. g. Q). 

F. One quasi proton configurations 

The lowest proton configurations are generated by oc- 
cupying the orbitals e, a and c in fig. ^ They are all of 
DAL type and < 80°. Accordingly, the configurations 
[e], [a] and [c] are interpreted as the AI — 1 bands with 
K'^ = 7/2", 7/2+ and 5/2+, respectively. The configu- 
rations [f], [b] and [d] are discarded. The configuration 
[g] has always z9o = 90°, as can be expected from fig. 
0. It is interpreted as the AI — 2 band (—,1/2), i. e. 



the favored signature sequence of the /ig/2 orbital. Fig. 
[T?] shows the calculated and the experimental Routhians 
in 7pTaio2- AH bands show the neutron AB crossing at 
UJ — 0.3 MeV . The TAG calculation for the configura- 
tion [g] gives too high energy and shows too early the 
neutron AB crossing. This is a well known problem of 
the hg/2 band which has been discussed in the literature. 
The discrepancies can partially be attributed to a larger 
deformation. Since these questions have been addressed 
before |6^,0 and are not at the focus of this paper we 
have not tried to improve the agreement by optimizing 
the deformation. 



G. One quasi proton one quasi neutron 
configurations 

The Routhians E'{d) for the four combinations of the 
quasi protons a and b emanating from Nilsson states 
[404] 7/2 with the quasi neutrons E and F emanating from 
[5 12] 5/ 2 are shown in fig. |l^. They are nearly degener- 
ate aX d — 90°. The configuration [aE] has its minimum 
at §0 = 35° and represents the AI — 1 band K'^ = 6^. 
The configuration [aF] has its minimum at 'do = 78° and 
represents the A J = 1 band K"^ = 1^. The other two 
configurations [bF] and [bE] continue the mirror images 
of [aE] and [aF]. They are discarded as spurious states. 
Both AI — 1 bands are seen in ^ ''^Lu. Fig. |l^ shows that 
separation and the slope of the bands is well reproduced 
by the TAG calculation. Another bundle of Routhians 
are four combinations of a and b with A and B, emanat- 
ing from the neutron states [633]7/2. The configuration 
[aA] represents the AI — 1 band ~ 7+ . The configu- 
ration [bB] has no minimum, only the kink at = 90°. It 
continues the mirror image of [aA] and is discarded. The 
configurations [aB] and [bA] have both their minimum 
at 'do — 90°. One is interpreted as a the AI = 1 band 
— 0+ and the other one is discarded. It does not 
matter which is chosen, because the two configurations 
differ only by their orientation {\aB >= e^'^^'^'^\bA >). 
To be definite we choose [aB]. 

When two quasi particles of the DAL type are com- 
bined into a low-_ftr and a high-X configuration , one must 
switch from the TAG to the PAG interpretation for both 
configurations simultaneously. As discussed in detail in 
ref. interpreting one configuration as PAG and the 
other one as TAG makes them highly nonorthogonal: The 
TAG solution is of the type ipKii^K2 whereas the PAG so- 
lution is of the type {^JRi iip-Ki){ipK2 ± V'-if2)/2- This 
has the following consequences for our example: 

i) One cannot take [aB] at its minimum at = 90°, 
because there it is of PAG type with good signature and 
thus nonorthogonal to the K — 7 configuration [aA]. It 
suffices to take the configuration at a somewhat smaller 
tilt angle, for example at = 85°, where [aB] has changed 
to the K = configuration . Energywise this makes 
barely a difference. However it is important for the cal- 
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culation of the B{M1) values, which according to ( |4l| , 
^) are zero for a PAC configuration . 

ii) Since the K"^ = configuration [aF] has a larger 'do 
than the 6~ combination [aE], its minimum approaches 
90° at a lower value of oj. However, one must keep the 
TAG interpretation as long as [aE] is tilted. As discussed 
for [aB], one may take i) = 85°. In fact, substantial Ml 
transitions are seen in this band, which are illustrated by 

Fig. |19| compares the experimental Routhians in 
Luio3 with the TAG calculations. The relative po- 
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sitions and slopes are well reproduced. Within the fre- 
quency range, all configurations are of TAG type. For 
the experimental K'^ = 0+ band [aB] the two signatures 
are separated (signature splitting). This is at variance 
with the calculations, which assign a A/ — 1 band (cf. 
preceding paragraph). The discrepancy is attributed to 
the residual interaction, which may lead to signature de- 
pendent correlations in the low-K bands, as well as to 
the zero point motion in d. 

It is noted that in yfTaioy the = 9~ and 0~ bands 
composed of the quasi neutron [624] 9/2 and the quasi 
proton [514]9/2 both do not show signature splitting [ |7^ ] 
and both have a substantial Ml transitions. This is con- 
sistent with the predictions by TAG. 

The orbitals H and G emanat e fro m the pseudo spin 
singlet [521]l/2 (cf. subsection |l^). The K"" = 4" 
configuration [aG] and the K'^ — 3^ configuration [aH] 
correspond to the parallel and anti parallel orientation 
of the pseudo spin with respect to the rotational axis. 
Accordingly, the distance between the 4^ and 3~ bands 
is equal to lo in the TAG calculation. The experimental 
distance somewhat deviates from this value, what can 
be seen as evidence for a pseudo spin dependence of the 
proton neutron interaction. 



H. Two quasi neutron excitations 

Fig. |l^ shows the Routhians E'{uj,-d) of the lowest 
positive parity configurations. Let us start with uj = 
0.2 MeV . For d > 60°, the first three configurations 
are [0], [AB] and [AG], which have the character of g- 
, s- and t- configurations, respectively. For both the s- 
and the t- configurations the ji components of the two 
quasi neutrons add up to a large value of Ji. In the 
case of the s-configuration the two js components are 
opposite in sign, resulting in J3 « 0, whereas for the t- 
configuration the two ja components add up to a large J3. 
The structure of the t- and s-bands was first discussed in 
where illustrations can be found. 
"For Lu = 0.2 MeV, the configurations [AB] and [AG] 
change order at 1} = 60°. For uj = 0.3 MeV, they mix 
strongly around z9 = 65°, interchanging their character. 
This refiects the quasi crossing between the orbital B and 
G, which can be seen in figs. p| and ^ For uj = 0.4 MeV, 
the crossing feature has disappeared. The orbitals B and 



G are now substantially different from what they were at 
low frequency. The labeling of the quasi neutron trajecto- 
ries suggested in figs. |o[ always assigns [AB] to the lower 
and [AG] to higher configuration . For oj — 0.2 MeV this 
results in an abrupt exchange of the structure when the 
configurations cross at 1? = 60°. [AB] takes the char- 
acter of a t-band for •& < 60° (which changes into the 
= 8+ state at = 0). Their sudden exchange be- 
comes a smooth transition at high frequency. 

The change of configuration [AB] from the s- to the t- 
structure reflects the strong response of the ^13/2 quasi 
neutrons to the inertial forces, which depend on the ori- 
entation of the rotational axis. In fig. ^ it is seen for 
the quasi neutron B as the rapid change of js from neg- 
ative to positive values near 60°. This is an example of 
the complex rotational response which cannot be guessed 
within the traditional PAG scheme. 

As seen in fig. ^ [AB] has its minimum at = 90°. 
It has the signature a = and is interpreted as the even 
spin s - band. It becomes yrast after the AB crossing 
at w = 0.3MeV. The configuration [AG] has also its 
minimum at 'i? = 90°. Since its signature is a = 1 it 
represents an odd spin band. The t-character of [AB] 
is only explored when additional quasi particles of DAL 
type change to smaller values. This will be discussed 
in subsection [III. For neutron numbers N > 106 the 



t-configuration is more favored, becoming a stable min- 
imum. The pure two quasi neutron t-band is seen for 
i80AAf.„„ iFnll and 
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Now we consider the combinations of the ^13/2 orbitals 
A, B, G, D with E and F, emanating from the Nilsson 
state [512]5/2. Fig. |2| shows the Routhians. In the 
lower bundle the quasi neutrons E and F are combined 
with A and B, emanating from [633]7/2. The configura- 
tion [AE] is the K'' = 6" band and [AF] the 1" band, 
both being A/ — 1 sequences. The configurations [BE] 
and [BE] are discarded. At w = OAAIeV the minimum 
of [AE] has moved to 80° and we switch to the PAG in- 
terpretation. Now, [AE] and [BF] have (tt, a) = (— ,1), 
i. e. they represent two odd spin bands, and [AF] and 
[BE] have (tt, a) = (— , 0), i.e. they represent two even 
spin bands. Fig. compares the calculation with the 
experimental bands in Y2'*Hfi02- Whereas the 6~ band is 
seen as a A/ = 1 sequence, as expected, the experimental 
1~ band shows a substantial signature splitting. The dis- 
crepancy is attributed to the residual interaction. For the 
low-i^ negative parity configurations the octupole corre- 
lations are important. Usually they are stronger for the 
(tt, a) = (-, 1) bands than for the (-, 0) bands ]0. This 
can explain why the experimental (—,1) sequence has a 
low energy relative to the experimental (—,0) but also 
relative to the TAG calculation. 

The interpretation of the bundle formed combining the 
quasi neutrons E and F with G and D is less straightfor- 
ward. The configuration [DE] has a TAG minimum for 
UJ < Q.25MeV. The component J3 w 2, i. e. it represents 
the 2^ band (cf. fig ^2|). Some small signature splitting 
is seen in experiment. 
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For Lu > O.SMeF we find d w 80° for [DE] and [CE]. 
Hence, the PAC interpretation is applied to the whole 
bundle. Accordingly, the TAG configuration [DE] splits 
into the PAC configurations [DE] and [DF], which have 
(7r,a) = (-,0) and (-,1). The TAG configuration [GE] 
splits into the PAG configurations [GE] and [GF], which 
have (—,1) and (—,0), giving rise to two A/ = 2 se- 
quences with odd and even spin, respectively. As seen 
in fig. there is a (—,0) band observed, which can be 
interpreted as [GF]. The odd spin band [GE] should be 
nearby. It is not given in the experimental level scheme, 
but ref. [^ reports two unplaced A/ = 2 sequences. 

The interpretation of the low frequency part of [GE] 
and [GF] has the same problems as discussed for the 
configuration [G] in sect 



[HE 



As seen in fig. gl| for 
Q.2MeV, [GE] is rather flat. It has a minimum 
at i9 = 0°, which is due to the K = 9/2 component in 
G. In addition it has a shoulder around t9 = 60° due 
the K — 5/2 component, which is hardly visible. How- 
ever the condition of uniform rotation (^) has solutions 
■& = 58° and 66° for w = 0.2 MeV and 0.25MeV, re- 
spectively. These points are included in fig. p^. The 
experimental (—,0) band, which we assign to [CF], is 
seen to substantial lower frequency than uj = 0.2, below 
which no TAG solution is found. This is another exam- 
ple for the limitations of the TAG approach, which treats 
the orientation angle as a static variable. The static 
approach is expected to work best if the function E{{}) 
is relatively symmetric around the minimum. Then, the 
fiuctuations of are also expected to be symmetric and 
the contribution linear in — i?o will average out. 

The combinations of the 113/2 quasi neutrons with the 
[514]7/2 orbitals will not be discussed, because the results 
of the calculations are similar to the ones obtained for the 
5/2~ orbitals and there is no experimental information 
about them. 

The combinations of the quasi neutron E with G and H, 
which emanate from the pseudo spin singlet [521] 1/2, are 
shown in fig. |2^. The situation is analogous to y^'^Luioa, 
where the same orbitals combine with the quasi protons a 
and b. The decoupled pseudo spin just adds or subtracts 
one half unit of angular momentum to the total angular 
momentum . This means, 



[EG] 



E', 



- Lu/2 + const, 
^\EH] = ^\E] + + const. 



[E] 



(72) 



The Routhians of [AG] and [AH] are related to the 
Routhian of [A] in a similar way (cf. figs. Qq and 23). 



The configurations [AG] and [EG] are observed as the 
K"^ = 4~ and 3+ bands. It would be interesting to 
observe the configurations [AH]3~ and [EH] 2+ with the 
pseudo spin being anti parallely oriented with respect to 
the rotational axis. Then one could investigate to what 
extend the the residual interaction depends on the pseudo 
spin orientation. 



I. Quasi neutron configurations at a large tilt 

The strong coupling estimate for the tilt angle 
■d = arccos(i^//) shows that the number of steps in / 
needed to change 1? from 0°( band head) to 90° increases 
with K. Hence, the high-iiT bands offer the possibility 
to explore the quasi particle spectrum at an angle sub- 
stantially different from 0° and 90°. The bands with low 
or moderate K, discussed in the preceding sections, only 
permit a sketchy view, because the rotational axis rapidly 
reorients from the 3- to the 1- axis. For i'i'4,i75jjj: ^^j^gj-g jg 
the family of bands built on the — 8^ proton config- 
uration [ae]. They permit us studying the quasi neutron 
configurations at a large tilt. These high- A' bands will 
be discussed in terms of the quasi neutron spectrum for 
1) — 45°. Although the tilt angle varies along the bands, 
it stays below 60° within the considered frequency range. 
Like the quasi particle diagrams for "i? ~ 90° (upper panel 
of fig. |l^ give a first orientation of the structure of the 
low-K bands, the diagram for -d = 45° (lower panel of 
fig. p^ ) permits qualitative interpretation of the quasi 
neutron configurations in the [ae] family. 

Comparing the two panels of fig. |l^ one notices a sub- 
stantial reordering of the trajectories. As a function of 
LU, they cross much more sharply for — 45° than for 
90°. This refiects the tilt of the rotational axis towards 
the symmetry axis, where all levels cross sharply. The 
sharp crossings have the consequence that the zero quasi 
neutron configuration [ae]8~ remains unperturbed up to 
the relatively high frequency of w « OAMeV, where the 
quasi neutron A+ interacts with the quasi neutron G'. 
This is confirmed by the full TAG calculation and the 
experiment shown in fig. The trajectory A+ crosses 
with a number of quasi neutron trajectories before it in- 
teracts with G'. This means that there are several two 
quasi neutron configurations below [ae]8~. The lowest 
of these, [aeAE], is seen in ^^"'Hf as the 14+ band. The 
lowest configuration of the opposite parity, [aeAB]16~, 
has not yet been observed. It is the quasi neutron t - 
configuration, which is shifted to low energy by the two 
DAL quasi protons a and e. Fig. ^ shows that with 
increasing oj it interchanges its character with [acAG]. 
This is the manifestation of the quasi crossing between 
B and G seen in fig. || at = 65°. The contrast to the 
quasi neutron spectrum at z9 « 90° is noted. There, the 
configuration [AB] interchanges gradually character with 
the vacuum [0]. For substantially smaller i), [AB] has a 
t - structure which is very different from [0] ( J3 « 8 and 
0, respectively) and, as a consequence, couples only very 
weakly to [0] . 

The 18"*" band [aeAEGI] is predicted too high relative 
to [aeAE] 14+. As discussed in section IIIJ, the discrep- 
ancy is probably due a substantial reduction of the pair 
field by blocking the two quasi neutrons G and I. 

The quasi crossing of A+ with G', seen oj ~ OAMeV 
in the lower panel of fig. |o|, causes the down bend of 
the zero quasi neutron configuration [ae]8~ seen at lOc = 
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0.38 MeV in fig. |^ It is the AB crossing at the tilt 
angle of i?o = 60°. It is delayed as compared to the yrast 
band where it appears at lJc — 0.30 MeV . Only part 
of the delay can be explained by geometry within the 
standard PAC scheme: The projection of the angular 
frequency on the 1 - axis is Wci — sin?? — ujc sin 60° = 
0.33 MeV , which is larger than loi = 0.30 MeV, where 
the AB crossing appears in PAC. 

Let us now discuss the [ae] family in the TV = 103 
system by means of the = 45° quasi neutron diagram 
(lower panel of The one quasi neutron configuration 
[aeA]23/2~ is lowest. It is observed at about the right 
energy in ^^^Hf as the 23/2~ band. Since E lies above 
A the configuration [aeE] is expected at higher energy. 
This confirmed by the full TAG calculation shown in fig. 
p5| . In the i} = 45° diagram the two quasi neutron 
excitations EB and EI have negative energy above cu = 
OAMeV. Fig. |2^ demonstrates that after minimizing 
with respect to d the configurations [aeABE]39/2+ and 
[aeAEI]35/2- he below [aeA]23/2-. This order of the 
bands is seen in ^^^Hf. 

Hence, the experimental high K spectra clearly reflect 
the modification of the quasi neutron spectrum with de- 
creasing z9. 

J. Unpaired neutron configurations 

In this subsection we discuss, how to interpret the high- 
K bands of the [ae] family in terms of configurations of 
the unpaired neutron system. Fig. ^ shows the single 
neutron Routhians a,t d = 45°. The labels of the states 
are chosen such that they are as close as possible to the 
quasi neutron neutron labels in fig. |l^. We consider 
N = 102. In order to keep the notation similar to the 
case of finite pairing, let us denote by [0] the yrast con- 
figuration for Lj < 0.1 MeV ( below the crossing between 
the levels A and E ). At z9 = 90° it has the character 
of the s - configuration and at d — 0° it is the ground 
state. The particle hole excitations are constructed rela- 
tive to this configuration [ae0]8~. Above the AE - cross- 
ing, the yrast configuration is [aeA^^E]14+. The next 
higher configuration is [aeA~^G~^EI]18+. They form the 
two lowest bands of the [ae] family. Both are shown in 
fig. 1^ as results of the full TAG calculation. The rel- 
ative position and slopes compare well with the experi- 
ment shown in the upper panel of fig. Glose by there 
is [aeA~^G^^BE]19~, which has not yet been observed. 
Above uj — 0.35 MeV one expects a number of config- 
urations with higher K, generated by lifting a neutron 
from the levels 9/2^ and G to I and B. One example is 
[ae A-iG-iEI]20+. 

The unpaired configurations [aeA^^E]14+ and 
[aeA~^G~^EI]18'*" are also shown in fig. ^for a compar- 
ison with the calculation at finite neutron pairing. There 
they are labeled in terms of the quasi particle notation 
as [aeAE]14+ and [aeAEGI]18+, respectively. The un- 



paired configuration [aeA~-'^E]18"'" lies significantly below 
the paired one. In this configuration four quasi neutrons 
are blocked and the description as unpaired neutron state 
is better (within the HFB scheme to which we restrict 
here). For [aeAE]14+, the paired calculation is favored 
at low and the unpaired at high uj. The zero pairing cal- 
culation compares better with the experiment. The zero 
pairing Routhians of both configurations lie too low as 
compared to the yrast band [0]. This may be a conse- 
quence of the dynamical pair correlations which are not 
taken into account. 

Fig. E8| shows the branching ratios for two bands 
[aeAE]14^and [aeAEGI]18+. In case of the = 14+ 
band, both the paired and the unpaired calculation give 
similar results. In case of the 18+ band, the unpaired cal- 
culation shows a similar increase at low frequency as the 
experiment, however it underestimates the experimental 
ratio. We have to underline here that the microscopic 
calculation of the magnetic transition probabilities by 
means of (^ ^) is expected to be less accurate than the 
popular strong coupling estimates (if applicable), which 
are based on quasi particle g - factors that are adjusted 
to the experiment. 

The unpaired configuration [ae0]8~ lies slightly below 
the paired configuration [ae]8~ after the latter has bend 
down. It has the character of the neutron s - configura- 
tion . That is the pairfield becomes small when the pair 
of ?i3/2 quasi neutrons is broken in the AB crossing. 

The unpaired configuration [aeA^^B]8^ lies above the 
two paired bands [acAG] and [aeAB], which interchange 
character. Thus the interpretation in terms of a quasi 
particle structure is favored. It will be interesting to see 
if the experimental bands shows the paired or unpaired 
pattern, which are markedly different. 

As seen in fig. ^6|, the lowest configuration in the 
N — 103 system is [ae A^^EI]35/2^, which is shown 
in fig. m as [aeAEI]35/2", A = 0. It lies below the 
paired configuration [aeAEI]. This is expected because 
three quasi neutrons are blocked, destroying the static 
pair field. The relatively high band head frequency in 
experiment is consistent with the zero pairing calcula- 
tion. In the calculation with finite pairing the band 
starts at a much lower frequency The configuration [ae 
A-iEF]23/2- is shown in fig. ^ as [aeA]23/2-, A = 0. 
It lies above the paired band [aeA]23/2^ for most of the 
frequency range. Finite pairing is favored for the one 
quasi neutron band. The zero pairing solution wins only 
at high frequency where a band crossing (down bend) 
is encountered. The band [aeAB I] 39/2+ corresponds in 
the unpaired scheme to [ae A"^BI]39/2+. This config- 
uration , shown as [aeABI]39/2+, A = in fig. |2^, lies 
above the finite pairing calculation. This somewhat sur- 
prising result (static pairing for the three quasi neutron 
configuration ) is understood from fig. In contrast to 
[ae A^^EI]35/2~, where the EF is blocked for pair scat- 
tering, for [ae A-iBI]39/2+ the pair BB' blocked. Since 
BB' is much further away from the Fermi - surface than 
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EF, blocking is much less effective. 

For higher frequency it becomes favorable to generate 
angular momentum by exciting quasi protons . As seen 
in fig. H, the next pair at large tilt is [eg]. In fact, the 
configuration [aecgAEGI]21^ appears at the yrast line. 
In the TAG calculation the paired four quasi proton con- 
figuration [aecg] has a lower energy than the unpaired. 
It is not unexpected that the TAG calculation gives a 
too high energy, because already the one quasi particle 
confi guration [g] is predicted too high by TAG (cf. sect. 



IV. RULES FOR TAG 

Let us summarize the experience gained in applying 
the TAG approach to the analysis of rotational spectra 
in axial well deformed nuclei. It is assumed that for = 
the rotational axis coincides with the symmetry axis. We 
shall use the form of rules. 

1. In order to construct the configurations use the 
particle - hole scheme in the case of zero pairing. 
For finite pairing use the quasi particle occupation 
scheme. The quasi particle levels appear in pairs 
of conjugate levels, one of which is occupied and 
the other is empty. Only for = 90° the conjugate 
levels have opposite signature. 

2. In order to check if a configuration corresponds to 
even or odd particle number trace it diabatically 
back to low uj, where gap between the quasi particle 
levels of positive and negative energy exists. 

3. In searching the equilibrium orientation -do try to 
stay within structurally the same configuration. 
Use diabatic tracing. 

4. If there are avoided crossings between the levels, 
diabatic tracing may end up in an unwanted con- 
figuration. Usually this shows up as an irregularity 
of the calculated quantities as functions of w. In 
such cases one has to resort to the quasi particle 
diagrams and manually reassign the desired config- 
uration . 

5. Reduce the chance of unwanted configuration 
changes by choosing the start angle of the dia- 
batic tracing close a to the expected equilibrium an- 
gles -do- Draw the quasi particle diagram e^{LO,'ds)- 
Never start diabatic tracing at = 90°, use 85°. 

6. The crossings between quasi particle configurations 
represent crossings between real bands, but the de- 
scription of the mixing region itself is incorrect. 

7. Rotational bands correspond to a function J{uj) 
which increases with uj. So long as = the band 
has not started yet. The band head lies at the fre- 
quency where a minimum at i?o 7^ appears. 



8. Solutions with do < 80° are of the TAG type. They 
represent A/ = 1 bands. The two signature part- 
ners (tt, a) and (tt, a + 1) are degenerate. 

9. Solutions with do > 80° are of the PAG type. They 
represent A/ = 2 bands of given (tt, a), i. e. I — 
a mod 2. The signature a is given by its value at 
d = 90°. 

10. When the tilt angle z?o(w) becomes larger than 80° 
the change from the TAG to the PAG interpretation 
results in unphysical jumps of the energy distance 
between signature partners and of other quantities. 
The experimental quantities show a gradual tran- 
sition between the two cases, which cannot be cal- 
culated by TAG. 

11. Since the number of AI — 2 bands of given (tt, a) 
must be the same in the PAG and TAG interpreta- 
tion, one half of the TAG configurations is spurious. 

12. Each configuration in the region < 90° has its 
mirror image E'{tt — — E'{{}) in the region i9 > 
90°. The diabatic continuation of the mirror image 
of an adopted configuration into the region i? < 90° 
is spurious and must be discarded. 

13. The spurious configurations have minima or kinks 
only at i9 = 90°. If such a configuration has another 
minimum for d < 80° this must be considered as 
physical. Then configuration has changed its char- 
acter with "iJ, being no longer spurious. 

14. For a strong tilt (0° < < 90°), the spurious 
configuration are usually high in energy and do not 
interfere. 

15. In the case of multi quasi particle configurations 
there are bundles of configurations emanating from 
■d — 90*^, each of which has its own do- Only when 
lowest do has reached 80° one must change from 
the TAG to the PAG interpretation simultaneously 
for the whole bundle. 

16. The intra band matrix elements of the electromag- 
netic operators can be calculated. PAG solutions 
provide their signature dependence. TAG solutions 
give only the average over both signatures. 

The rules can be applied to triaxial nuclei with few 
obvious modifications. Rule 5 must be complemented by: 
Do not use dg ~ 0°, start at 5°. Only the first, general 
statement of rule 7 remains valid. The second which 
assumes that i9 = is a symmetry axis does not apply. 
The low frequency behavior of triaxial nuclei has not yet 
been studied, except the investigation of a model case in 
[|75| . PAG solutions are possible for all three principal 
axes and TAG solutions are possible in all three planes 
spanned by the principal axes. Accordingly one must 
extend the search. The simplest way is letting 7 vary 
from 120° to -60° (cf. tab. 1). If the rotational axis does 
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not lie in one of the three planes spanned by the principal 
axes the rotation becomes chiral. The consequences of 
this 3D geometry are discussed in and |7|] . 

V. CONCLUSIONS 

The semiclassic concept of uniform rotation about an 
axis that is tilted with respect to the principal axes of the 
deformed density distribution leads to a mean field theory 
which describes energies and intra band electro-magnetic 
matrix elements of A/ — 1 bands in a quantitative way. 
The orientation of the rotational axis turns out to be as 
good a collective degree of freedom as the familiar shape 
degrees of freedom are. The tendency of high spin par- 
ticles to align with the rotational axis, which in general 
does not have the direction of one of the principal axes 
of the deformed mean field, is a concept that permits 
to explain many features of high-ii' bands from a new 
perspective. 

The tilted solutions do not have the familiar C2 sym- 
metry, which appears when rotational axis coincides with 
one of the principal axes. The lower symmetry results in 
the loss of the signature quantum number, which mani- 
fests itself by the appearance of one AI = 1 band instead 
two separate AI = 2 sequences of opposite signature. 
The transverse magnetic dipole moment, which deter- 
mines the rate of magnetic dipole transitions, plays the 
role of an order parameter. For tilted solutions it has a 
finite expectation value, which may become quite large, 
because it is the sum of contributions of several quasi 
particles . For rotation about a principal axis the expec- 
tation value of the transverse magnetic dipole moment is 
zero. The magnetic transition probability is given by a 
matrix element between two different quasi particle con- 
figurations , which is of single particle order. 

Being a mean field approach, tilted axis cranking is 
not capable of describing the transition from a tilted to 
a principal axis solution in a correct way, because this 
involves the transition from a broken to a restored dis- 
crete symmetry. The signature dependence of the en- 
ergy and other quantities appears in a sudden, unphysi- 
cal way when switching from the broken symmetry to the 
conserved symmetry interpretation. Still, one can guess 
from the calculations at which rotational frequency the 
signature effects are expected and how strong they should 
be. 

The breaking of the C2 symmetry leads to the appear- 
ance of spurious states. An elimination method is sug- 
gested. After applying it the calculated sequence of the 
first bands above the yrast line which agrees with the ob- 
served one for the studied examples. No spurious states 
remain in the near yrast region. 

It is the strength of the tilted axis cranking approach 
that treating many excited quasi particles is no more 
complicated than treating few. In order to demonstrate 
the application of the method we studied configurations 



with up to four excited quasi protons and four excited 
quasi neutrons in the nuclides with N — 102, 103 and 
Z ~ 71,72,73. The calculated energies and branching 
ratios agree with the experimental values within an ac- 
curacy that is typical for microscopic mean field calcula- 
tions. In particular, it is found that the order and struc- 
ture of the high-X bands can be qualitatively understood 
in terms configurations constructed from quasi particle 
levels, which are calculated as functions of the rotational 
frequency w at a fixed tilt angle of about 45°. 

The regime of quenched static pairing is encountered in 
the multi quasi particle bands of high seniority. Since the 
change of the pair field is not at the focus of this paper, it 
was treated in a rough way by comparing the paired with 
unpaired solutions and choosing the one with the lower 
energy. This schematic approach turned out quite practi- 
cable for a first analysis of the high-X band structure. A 
refined description of pairing within the tilted cranking 
model, which includes dynamical pair correlations, will 
be given elsewhere Q . 

The study of high- if bands at the largest frequencies 
attainable is a interesting problem of nuclear physics. 
More systematic studies than the present are expected to 
reveal the response of the single particle levels to the tilt 
of the rotational axis. A particularly interesting question 
is how the K quantum number is eroded with increasing 
rotational frequency. Tilted axis cranking is the proper 
mean field theory to address this question. It is also the 
appropriate starting point for theories that go beyond 
the mean field, like RPA. 
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FIG. 2. Quasi proton energies for e = 0.258. £4 = 0.034, A = 0.75A4eV . The chemical potential A is adjusted to have 
< Z 72 at w = 0. Full lines: positive parity, dashed-dotted lines: negative parity. For ■& = 90° the signature is indicated 
by ± standing for a = ±1/2. 
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FIG. 3. Quasi neutron energies for e = 0.258, £4 = 0.034, A = Q.GQMeV . The chemical potential A is adjusted to have 
< N >— 102 at = 0. Full lines: positive parity, dashed-dotted lines: negative parity. For 1? = 90° the signature is indicated 
by ± standing for a = ±1/2. 
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FIG. 4. Routhians of the configuration [0] for different values of the pair gaps (upper panel) and of the chemical potential 
A (lower panel). The quantity E' + X^Z + X^N is displayed. The pairs of numbers give (Aj,, A„) in MeV. The values 
Ap = 43.95 MeV and A„ = 49.0 MeV, corresponding to < Z >= 72 and < N >= 102 for a; = 0, are used in the upper 
panel. In the lower panel Ap = 43.95 MeV and A„ is explicitly indicated in MeV. The arrows show at which u! the condition 
< AT >= 102 is satisfied. 
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FIG. 5. Single proton energies at intermediate frequency for e = 0.24, £4 = 0. Full lines: positive parity, dashed-dotted 
lines: negative parity. For -ff = 90° the signature is indicated by ± standing for a = ±1/2. On the right hand side the Nilsson 
labels are given, which are relevant for to = 0. 
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FIG. 6. Single proton energies at high frequency for e = 0.24, £4 = 0. Full lines: positive parity, dashed-dotted hnes: 
negative parity. 
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FIG. 7. Single neutron energies at intermediate frequency for e = 0.24, £4 = 0. Full lines: positive parity, dashed-dotted 
lines: negative parity. For § = 90° the signature is indicated by ± standing for a = ±1/2. 
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FIG. 8. Single neutron energies at high frequency for e = 0.24, £4 = 0. Full lines: positive parity, dashed-dotted: lines 
negative parity. 
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FIG. 9. Angular momentum components of different quasi neutrons as function of the tilt angle i? at = 0.3 MeV. They 
are labeled by A, B, ... as in fig. H. The dotted lines show the values ja = K and ji = ji{uj = 0.3 MeV s'm^,^ = 90°), 
which are used in the standard CSM treatment of bands with finite K. They are labeled by AP, BP, ... The parameters of the 
calculation are same as in fig. H and pi 



31 




32 




FIG. 11. Quasi neutron energies for e = 0.258, £4 = 0.034, A = OmMeV and uj = 0.4 MeV. The chemical potential A is 
adjusted to have < N >= 102 at w = 0. Pull lines: positive parity, dashed-dotted lines: negative parity. 
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FIG. 12. Total Routhian as function of the tilt angle for the one quasi neutron configurations in 72 Hfios. 
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FIG. 13. Total Routhian E'(uj,'&)/uj as function of the tilt angle for the one quasi neutron configuration [E] in 72^Hfio3. 
The gridlines correspond to the steps Ai? = 5° and Aw = 0.02 MeV, starting with w = 0.02 MeV. A constant is added such 
thati;'(a;,i? = 90°)/a; = 0. 



35 



80" 



60° 



40° 



2(f 



0° 








0.3 



0)[MeV] 



FIG. 14. Equilibrium tilt angle ■&„ as function of the rotational frequency uj for the one quasi neutron configurations [A] and 
[E]. The full lines display the TAG result. The dashed lines show the strong coupling estimate, where the moment of inertia 
J' = 35 MeV~^ is used, which is the value of J/uj for the configuration [E]. 
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FIG. 15. Total Routhians as function of the tilt angle for the lowest positive parity quasi neutron configurations in 72 Hfio2. 
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FIG. 16. Total Routhians of the lowest one quasi neutron configurations in 7|^Hfio3. 
The following legend applies to all figures of this type and contains more information than needed for this figure. 
The line convention for the parity and signature (vr, a) is as follows. 
Odd A: full (+, 1/2), dotted (+, -1/2), dashed dotted (-, 1/2), and dashed (-, -1/2). 
Even A: full (-f, 0), dotted (-(-, 1), dashed dotted (-, 0), and dashed (-, 1). 

Upper panel: Experimental Routhians from the data in |]6^ , p6| . If the bands do not show signature splitting the Routhians are 
calculated by means of eqs. ( |5^ - |5^ ) if there is a signature splitting eqs. ( ^2[ - |53| ) are used. The configuration is given as [X]Ji''^ 
where X stands for the letter code explained in the text and K"^ gives the parity vr and the K value of the band head. 
Lower panel: Results of the TAG calculations. For TAG solutions (i9o < 80°) the two lines representing the signatures are on 
top of each other. If they are separated the solution is of PAG type (i9o > 80°). The last TAG point and the first PAG point 
are connected by a straight line. The grid is Aoj = 0.05 MeV . The band head is only determined with the accuracy of Acj. 
For several bands the frequency of the first transition ( J(t<-') = K ^\) is indicated by a fat dot. 
A rigid rotor reference with JT" = 50 MeV~^ is subtracted from all Routhians. 
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FIG. 20. Branching ratios of the one quasi proton one quasi neutron configuration [aF]j^- in yi^Luioa. 
Experimental values from [^. Full line: TAG calculations. 
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FIG. 23. Total Routhians of the two quasi neutron configurations containing the pseudo spin singlet quasi neutrons G and 
in 72*Hfio4. Upper panel: Experimental Routhians from the data in Lower panel: TAG calculations. Gf. caption of fig. 
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FIG. 24. Total Routhians of the lowest quasi neutron configurations in 72 Hfio2 built on the K'' = 8" two quasi proton 
excitation. 

Upper panel: Experimental Routhians from the data in ^]. Lower panel; TAG calculations. The labels A and denote 
A„ = 0.69 MeV and 0, respectively. Cf. caption of fig. 
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FIG. 25. Total Routhians of the lowest quasi neutron configurations in 72^Hfio3 built on the — 8~ two quasi proton 
excitation. 

Upper panel: Experimental Routhians from the data in ^]. Lower panel; TAG calculations. The labels A and denote 
A„ = 0.69 MeV and 0, respectively. Gf. caption of fig. 
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FIG. 26. Single neutron energies for i9 = 45°, e = 0.258, £4 = 0.034. 
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FIG. 27. Total Routhians of the lowest neutron configurations in ^rHfio2 built on the — 8 two quasi proton excitation. 
The TAG calculations assume zero neutron pairing. Gf. caption of fig. u]l 
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FIG. 28. Branching ratios of the configuration [aeAE]]^4+ and [aeAEGI]ig+ in 72*Hfio2. 
Experimental values from pa. Full lines: TAG calculations with neutron pairing dashed lines with zero neutron pairing. 
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TABLE I. Association of the principal axes of the triaxial potential with the axes 1 and 3. The shapes are labeled by p 
(prolate), o (oblate), and t (triaxial). The axes are labeled by 1 (long), i (intermediate), and s (short). The star indicates a 
symmetry axis. 
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